Mott Scattering of polarized electrons in a circularly polarized laser field 
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We present a study of Mott scattering of polarized electrons in the presence of a laser field with 
circular polarization using the helicity formalism and the introduction of the well known concept 
of non flip differential cross section as well as that of flip differential cross section. The results we 
have obtained in the presence of a laser field are coherent with those obtained in the absence of a 
laser field. We have studied the relativistic regime as well as the non relativistic regime that are 
precisely defined in the text. Two important consistency checks have been carried out successfully. 
The first one is that the sum of both differential cross sections (one with spin up, the other with spin 
down) always gives the unpolarized differential cross section. The second one is that the relativistic 
unpolarized differential cross section converges to the non relativistic differential cross section in the 
limit of small velocities. Moreover, the results obtained using the sofware Reduce [3| gave rise to 
non contracted coefficients that have been dealt with using the geometry chosen. 
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I. INTRODUCTION 

This work deals simultaneously with three important 
topics in Atomic Physics, namely the spin of a Dirac or 
a Dirac- Volkov particle, the concept of a spin polarized 
relativistic particle and finally how the laser field affects 
the process of scattering of such particles. We first give 
a brief overview of the laser. The word laser is by now 
well know to the lay-man because of its applications in 
the field of medicine ( surgical and diagnostic procedures 
), telecommunications ( fiber-optic telephone links, com- 
pact disk information storage, etc ) and technology ( laser 
drilling of materials, geodesies measurements, newspaper 
printing ). Lasers come in different shapes, sizes and 
prices and under different names as ruby laser ( the first 
laser to operate jlj ), helium- neon-argon, semi-conductor 
laser and others. It is not possible to explain what a 
laser is and how it works in a few words. At this stage, 
it suffices to say that a laser is a light source with pe- 
culiar characteristics, drastically different from those of 
conventional sources such as a candle or a light bulb and 
is radiated in a single direction ( not in all directions as 
in a light bulb ) which enables any person working in the 
field to collect it in a lens and focus it in a region of very 
small dimensions. The spectral purity of this process 
and the directionality of the laser light dramatically im- 
proves the efficiency of this procedure, making it possible 
to concentrate a sizeable amount of power in a small re- 
gion for different operations ( like the melting or cutting 
of metals.) The laser is basically used as a very powerful 
light bulb. However they are others ( like optical commu- 
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nications ), in which its most important characteristics 
are the spectral band- width and angular aperture of the 
emitted beam. To understand them, one needs to con- 
sider what light is and how it is emitted, which in turn 
depends on the emitter, the atom. This task requires 
the introduction of some basic concepts of quantum me- 
chanics. In this work, however, it is possible to describe 
the laser field classically since the numbers of photons in- 
volved is high 0. As for the spin of a particle, according 
to Landau and Lifshitz this property of elementary 
particles is peculiar to quantum theory and therefore has 
in principle no classical interpretation. In particular, it 
would be meaningless to imagine the "intrinsic" angular 
momentum of an elementary particle as being the result 
of its rotation about its own " axis" . It is a fascinating 
and mysterious complication. Moreover, its practical ef- 
fects prevail over the whole of Quantum Physics. The 
existence of spin and the statistics associated with it is 
the most subtle and ingenious design of Physics, without 
which the whole universe will collapse. Spin occupies a 
unique position since a wide range of Physics with differ- 
ent degrees of difficulty is needed for its understanding. 
The fact is that most textbooks ( even the most advanced 
ones ) devote at most one chapter to this subject and give 
only a utilitarian description. The theory of special rela- 
tivity, which is essential for the understanding of spin and 
statistics is often forgotten except in deriving the Dirac 
equation Q. The relation between spin and statistics is 
hard to understand. Feynman Q wrote : " it appears to 
be one of the few places in Physics which can be stated 
very simply, but for which no one has found a simple and 
easy explanation. The explanation is down deep in Rel- 
ativistic Quantum Mechanics 3. This probably means 
that we do not have a complete understanding of the 
fundamental principles involved " . Now, let us turn to 
the physics of spin polarized free electrons. Our aim is 
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to oppose the current trend in science that sees special- 
ists writing primarily for like-minded specialists. It is 
for example necessary to give a simple introduction of 
Mott Scattering one of the most important techniques in 
polarized electron studies. Only in recent years, has it 
been found possible to produce electron-beams in which 
the spins of the particles have a preferential direction 
0. There are many reasons for the interest in polarized 
electrons. One essential reason is that in physical investi- 
gations, one endeavors to define as exactly as possible the 
initial and/or the final states of the systems being consid- 
ered and also to have electrons in the form of a well de- 
fined beam, that is a beam in which the directions of the 
momentum of the individual electrons are as uniform as 
possible. Polarization effects in electron scattering were 
assumed to be significant only if the electron velocity was 
comparable to the velocity of the light. It was not un- 
til the 1960s that large polarization effects in low-energy 
electron scattering were ascertained [2]. It is important, 
before presenting our investigation about laser-assisted 
Mott scattering of polarized electrons to sketch the prin- 
cipal steps of our treatment. For purpose of clarity and 
simplicity, we begin by the most basic results of Mott 
Scattering of polarized electrons in the absence of a laser 
field. Then, in the presence of a circularly polarized laser 
field, we give a detailed account of the formalism used 
and we compare the results with those obtained in the 
absence of a laser field. The organization of this work 
is as follows : in section 2, we present the scattering of 
polarized electrons by the Coulomb field of a heavy static 
charge in the absence of a laser field and we present the 
concept of a polarized differential cross section ( in brief 
DCS ). We also present the helicity flip DCS as well as 
the helicity non flip DCS. Then, we define the degree of 
polarization of the scattered electrons. At this stage, it 
is of a paramount importance to remark that in exper- 
iments, the degree of polarization of the scattered par- 
ticles is measured '?|. In section 3, we present the laser 
assisted Mott scattering of polarized electrons in a laser 
field with circular polarization. After section 4, devoted 
to the discussion of the results obtained, we end by a 
brief conclusion in section 5. Throughout this work, we 
use atomic units h = m = e = 1 and work with the met- 
ric tensor g^'^ — g^^^ = diag(l, —1, —1, —1). The angular 
frequency we have chosen is = 0.043 (a.u)= 1.17 eV. 
It corresponds to the lasing transition of the Nd laser at 
a wavelength of 1064 nm. 



element corresponding to this process. 
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describes the incident electron normalized to the volume 
V and 



1 



(3) 



describes the scattered electron with the same volume 
normalization. The Coulomb potential is gen- 

erated by a static heavy nucleus of charge — Z, that is 



<o„z(^) = (-0'0) 



(4) 



Therefore, we have for 4-coui{x) : 

4coulix) = 7°Ao = - 



The transition matrix element is then given by : 
' c J J ^/2EjV |x| VWV^ 

(5) 

The integral over xq gives : 

r' + oo p-\-oo 



dx°e-'^P^o-Pfo)-^" = 

= c.2ndiEi - Ef) (6) 
while the integral over x gives : 

J (P» - Pf? |q|' 

where we have introduced the momentum transfer q = 
Pi — Pf- Using the standard procedures of QED we 
get for the unpolarized DCS. 



da Z2 1 1 



dQf c4 |q|4 2 



^ |m(p/,s/)7°u(p„s,)P (8) 



II. MOTT SCATTERING OF POLARIZED 
ELECTRONS IN THE ABSENCE OF A LASER 
FIELD. 



We begin by one of the simplest process of QED, 
namely the process of Mott Scattering of electrons in the 
lowest order of perturbation theory. We follow the usual 
steps of calculations and we give the transition matrix 



evaluated for Ei = Ef = E. As a side-result this implies 
in turn |pi| = |p/| = |p|. InEq.® we have averaged over 
the initial spin polarizations Si and summed over the final 
spin polarizations s f . We obtain the simple result : 



:2\r,EiEf 



(9) 
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and with Ei = Ej 
end up with : 

\ l"(P/'«/)7°"(Pi,s»)P 

Si,Sf 



E as well as |pi| = |p/| = |p|, we Therefore, the definition of a 4-spin vector follows : 



2c'[^ + |p|'cos(0)+c2] 



where d is the scattering angle and {pi.pf) = E'^/c 



(10) 



Pf.pi = E-^/c^ — \p\^cos{9). Using the relation 
c2|p|2 = H) and |q| = |p, - p/l = 2|p|sin(^/2), the 
final expression for the unpolarized DCS is : 

da 1 1- (3^sm'^{0/2) 



dfl 



f Mott 



4c2/32|p|2 sin4(0/2) 
= ^ (l-/3^sin2(0/2)) (11) 



/ Ruth 



where 



da 



d^fRuth 4c2/32|p|2sin4(0/2) 



(12) 



is the Rutherford unpolarized DCS in the limit of small 
velocities {(3 0). The calculations carried out up to 
now assumed that the initial electron spin is not ob- 
served. So, we focus now on the scattering of polarized 
Dirac particles. For that purpose, we need some formal- 
ism to describe the spin polarization Free electrons 
with 4-momentum p and spin s are described by the free 
spinors u{p,s). The spin 4-vector is defined by : 



c c 



(13) 



where p = p/|p| is the unit vector defining the direction 
of the 3-vector p. The 4-vector s'^ is a Lorentz vector in 
a frame in which the particle moves with momentum p. 
The following properties can be easily checked : 
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(14) 



This is the normalization relation for the 4-vector s^. 
One also have : 



A 



E 



5<: = -(ipi,-p) 



(19) 



We calculate now the polarized DCS for Mott Scattering 
of an electron with a well defined momentum pi and a 
well defined spin Si. If also the final spin s/ is measured, 
the polarized DCS reads : 



da _ Z'^ \u{pf,sf)'y°u{p^,s^)\''' 



f 



(20) 



In Eq.([20|). one has to be careful since there are no sum- 
mation on spin ( either initial or final ) polarizations. 
We introduce the two operators : 



Sa,(s/) =i(l + A/75^/) 
and we obtain for the polarized DCS : 



(21) 
(22) 



da 
d% 



(A., A/) 



Z^l 

c4 4 



Tr{j°{l + X,-f5^i){^^c + c^) 
7°(l + A/75^/)(;5/c + c2)} (23) 



using the relations (pi-Sf) = J-|p|(l - cos(6')) = ipf-Si), 
(si.Sf) = (IpP — cos(0)/c2) /c^ as well as (pi-Pf) = 

^-|p|'cos(0) 
we obtain : 

J^(A., A/) = |^2{^' + 

+X,XfiE^ cos2(6'/2) - sin2(0/2))| (24) 



At this stage, let us note that the Xi and A/ have the 
following properties : 



A? = A^ = 1 



(25) 



(p.s) 



-p^.s^^Q 



(15) 



It is the orthogonality relation between p and s. Next, 
we introduce the spin projection operator : 



A,;.Af = 1 



Xi=Xf = 1 
X,^Xf = -1 



(26) 



1 



S(s) = -(l+75^) 



(16) 



with 75 = ij'^j^^'^j^ 
simple property : 

E(s) u{p, +s) — u{p, +s) 



-170717273. This operator has a 



t{s)uip,-s)^0 (17) 



This formalism can be applied to helicity states where 
the direction of the spin points along the direction of the 
momentum 3-vector p. 



S,=X^ 



A = ±1 



(18) 



A,;. A 



X,=-Xf = 1 
A, = -A/ =-1 



(27) 



So, if during the scattering process. A/ = —Xi which 
means that a helicity flip occurred, the flip polarized DCS 
is : 



da 

dn} 



flip 



4Z' 



^sin2(6'/2) 



da 
dih 



Ruth 



E^ 



(28) 
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The case where there is no hehcity flip corresponds to 
Aj = Xf, so that XiXf = 1 and the hehcity non flip po- 
larized DCS is given by : 



da 



nonflip 



da 

dn^ 



(29) 



Ruth 



Of course, the sum of the hehcity flip polarized DCS and 
the hehcity non flip polarized DCS must give the unpo- 
larized DCS. 



da 



da 



dVtf 

Noting that E"^ 



Mott dQf 



da 



nonflip 



& + or 



dVtf 



(30) 

flip 

pI^c^ we have 



sin2(^/2) + 008^(^/2) = 1-^2 sin2(e/2) (31) 

E 

and the flnal results is : 



da 



dflf 



Mott 



da 



[l-/32sin2(0/2)] (32) 



Ruth 



This result is of a paramount importance. 

For every process, the sum of the flip polarized DCS and 
the non flip polarized DCS always gives the unpolarized 
DCS. In fact, this is theoretically exact but, as we shall 
see very soon (section 3), when it comes to numerical sim- 
ulations and where infinite summations over the number 
of laser photons arc involved, due to a unavoidable lack 
of precision, the theoretical exact result becomes a nu- 
merical approximate result. Another important result is 
the degree of polarization deflned by : 



P = 



da 

da 



'-iX, = Xf = l)-^{\ = -Xf = l) 



A/ 



-Xf = 1) 



(33) 



For our process, P reads 



P = 



cos2(6'/2) - |^sin2(6'/2) 



cos2(6'/2) 



1-2- 



^ sin' (0/2) 
sin' {9/2) 



E'' cos2 (61/2) + c4sin^ (61/2) 



and with E = ^c', one has : 



P = 1-2- 



sin2(6'/2) 



7^ cos^ (6'/2) +sin^(0/2) 



(34) 



(35) 



We notice that for the process of Mott scattering in the 
absence of the laser fleld, only two parameters are rele- 

'1-/3' and the 



vant : the relativistic parameter 7=1/ 
scattering angle Moreover, one can see from Fig. (1) 
that when the collision approaches the relativistic do- 
main, the degree of polarization becomes less strongly 
dependent on the relativistic parameter 7 and becomes 
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FIG. 1: Degree of polarization as function of the scattering 
angle 6if, for different values of the relativistic parameter 7. 



close to a constant value P ~ 1. The value P = 1 cor- 
responds to a collision when no spin flip occurs whereas 
the value P — —1 corresponds to the reverse situation 
(there is a spin flip). In Fig. (2), we plot P as a function 
of the final electron momentum 6f and a different pic- 
ture emerges. Since the scattering angle is not involved, 
the minimum value of P shifts to a value greater than 
P = — 1 and this trend is more and more noticeable for 
increasing values of 7. 



III. MOTT SCATTERING OF POLARIZED 

ELECTRONS IN THE PRESENCE OF A 
CIRCULARLY POLARIZED LASER FIELD. 

The 4-vector potential A^^ chosen is such that : 

A" = cos{4>) + sin(^) (36) 



with (ai.ai) = -\a\' = -\A\', (02. a2) = -\a\' = -\A\' 
and the polarization 4-vectors Oj* and ai^ are such that : 



< = |a|e^^ = |a|(0, 1,0,0) 



and 



ae- 



|a|(0, 0,1,0) 
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FIG. 2: Degree of polarization as function of the angle 6f, 
for different values of the relativistic parameter 7. 



In Eq.®, = (k.x) = fc^.x^ with fc^ = (fcO, 0, 0, k^) = 
Ci;/c(l, 0, 0, 1). satisfies the Lorentz condition : 



kf^.A'^ = 



(37) 



which imphes (fc.ai) — {k.a2) — i.e, we choose k along 
the direction of the Oz axes. 

Now, in presence of a laser field, the Dirac wave functions 
describing the incident and scattered electrons have to be 
replaced by the Dirac- Volkov wave function Q which are 
given by : 



%[x)^R{q) 



u(j>, s) 



where 



Riq) = Rip) = 1 



= 1 



(38) 



(39) 



2c{k.q) 2c{k.p) 
since the 4- vector that now replaces p^^ is such that : 



= pf" 



a' 



2c2(fc.p) 

The time component of is such that : 



Q Q E a llJ 1 a^uj 
c c 2c^{k.p) c c 2c'^{k.p) 



(40) 



) (41) 



We shall see that this 4- vector is the new canonical im- 
pulsion 4-vector that the electron feels inside the laser 
field. The 4-vector q^ is called the quasi-impulsion and 
has the following property : 



a^k 
2c\k.py^ 



[P- 



a^k 



2c^{k.p)' 



P,.p- - ^ 



Remark also that we can write : 

R{q) = [1 + c{qm] = [1 



(42) 



(43) 



with c{q) = c{p) = l/[2c(fc.g)] — l/[2c(fc.p)], since 
(k.q) = (k.p) and k^ = kf'.k^ = 0. The function S{q,s) 
appearing in Eg. ([38]) is such that : 



Siq, s) = -q.x sin(0) + cos(^) (44) 

One has for the transition matrix element : 
i 



S 



d xxl:g^{x)4couiix)ipq^ix) (45) 



One must be very cautions not to confuse the Coulomb 4- 
vector potential due to the charge of the nucleus with the 
4-vector potential of the circularly polarized laser field. 
As before one has 4-coui(x) — — Z7°/|x| and therefore : 



Z f ^ 



with 



and 



'ipq^ix) = R{pi) 



u{pi,s. 



^2QV 



we transform giS{qi,x)-iS{qf,x}^ jj- introduce 



with 



and 



ai 



1 f {ai.pi) _ jai.pf) 
(k.pi) (k.pf) 



^ ^ 1 / {a2-Pi) _ (Q2-P/) 
^ c V (k.pi) (k.pf) 



(46) 



(47) 



(48) 



(49) 



(50) 



(51) 



we get : 

S{qi,x) ~ S{qf,x) = ~{qi - qf).x - zsin(0 - (/jq) (52) 
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where : 



with 



(53) 



cos(0o) = — or sin((/)o) — — 



Now, let us look closely at the quantities in Eq. (|46l 
First, we have : 



In particular, we remark that : 

i?(p/)7"i?,(p,) = {1 + c(p/)[{^i^cos(/) + (i2^sin(/)]} 
X7°{1 + c{pi)[^^i coscj) + sincj)]} 
= {-f° + c{pf)[4i^j° cos0 + ^2|^7° sincj)]} 

X {I + c{pi)[^^i cos (j) + li^2 sin (j)]} (55) 

Before invoking the well known relations involving ordi- 
nary Bessel functions, we have first to transform Eq. (|55p 
in a more compact form. Eq. (|55p contains nine terms, 
and after some algebraic calculations the final result for 
i?((Z/)7°i?(g.) = i?(p/)7°i?(p») is : 

R{qfh"R{q,) = 7" - 2^^c(K)c(p/)a'|^ 

+{c{pih'^Hi + c{pf)4ih"} cos <j) 

+{c{Pth"H2 + c{pf)42h"}sin(t) 

= Co + Ci cosict,) + C2 sin((/)) (56) 



with 




^o_2^c{pMPf)a^^ 
c 

c{pi)-f°l^4i+c{pf)4il^j° 

c{p,)-i°U2+c{p})i2H 



iz sin{0— 0o) 



E 



Bn{z) 
Bln(z) 
B2n{z) 



-in(pi 



(57) 



Bn{z) = J„(z)e™"*« 

Bin{z) = (J„+i(z)e'("+i)^« + J„_i(z)e'("~i)'^'')/2(58) 
B2n{z) = (J„+i(z)e*("+i)^« - J„_l(z)e*("-l)*°)/2^ 



'4^qf{xh°'4'q,{^) ^ Mpf,Sf)R{pf)^°R{pi)u{pt,St) (54) Therefore: 



+ CXD 



R{pf)fR{p,) = [CoS„(z) + C2Bi„(z) 

n— — oo 

+C3S2„(^)]e-"'^ 



+ 00 



-in(f} 



(59) 



The transition matrix element becomes : 



iZ2TrS{Qi — Q f + nuj] 



47r u{pf,Sf) ^ u{p.„Si) 
|q|2 ^/2QjV "ySQlF 



(60) 

Using the standard procedures of QED, one obtains for 
the polarized DCS, we have : 



da 



+2? du^ 

1. — — ryn 



(61) 



where : 



da(") \qf\Z^ 1 2 
- — -X-74 l"W>s/)rnw(p„sO| (62) 



|q| 



evaluated for Q/ = + nw. 
The squared matrix element is : 



\u{pf,Sf)T,^u{p^, Si)p = ut(pi, Sj)rjj7°u(p/, Sf)u{pf,Sf)Tnu{pi, s,) 

= u{pi,s^)'y°Tl'y°u{pf,Sf)u{pf,Sf)Tnu{pi,Si) (63) 

If the operators T.Xi{si) and E>^(s/) are introduced, we obtain : 

\uiPf,sj)Tu{p., s.)f = Tr (r„ii±^:^(Ac + c^)Tjl±^pllli^,c + A (64) 



where : IV. RESULTS AND DISCUSSIONS. 

r„ = 7°rt7° 

_ j^o _ 2—c{p )c{p )a'^^]B* Before presenting the results and their physical interpre- 

c * " tation, we would like to emphasize the fact that the RE- 

+ [c(p,)!^i^7° + c(p/)7"^^i]Sr„ 
+ [c{p,)Ml' + 4Pfh°U2]B;^ (65) 
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DUCE code we have written gave very long analytical 
expressions which were difficult to incorporate in the cor- 
responding FORTRAN code we wrote to extract figures 
and tables. From now on, the pair of indices n, / will 
stand for " non flip " and the index / alone will stand 
for " flip ". For the non flip DCS, let us simply write 



da 

dh 



n.f 



E 



n— — oo 



(66) 



nj 



where 
da^' 



dflf 



(nJ) 



+B{X, = A/ = l)[J,^+i(z) + J^i(z)] (67) 
+C(A.j = Xf ^ l)J„+i(z)J„_i(z) 

+V{\, =Xf = l)J„(z)[J„+i(z) + Jn-l{z)]} 

The four coefficients A{Xi,Xf), B{Xi,Xf), C{Xi,Xf) and 
'D(Ai,A/) are very long to write down so we prefer to 
focus on their global contents instead of giving tedious 
details and explanations. First, we noticed that in both 
coefficients Anj = A{Xi = A/ = 1) and Cnj = C{Xi = 
Xf — 1), there is no occurrence of the various com- 
pletely anti-symmetric tensors £0/3^5 (where the indices 
are Lorentz ones and take integer values from to 3). 
This clearly means that these tensors were totally con- 
tracted and we ended up with two tractable coefficients 
that where very easy to incorporate into the main FOR- 
TRAN code. Let us remind the reader that we used 
throughout this work, the convention : 



£0123 



1 



(68) 



meaning that £ap^& — 1 for an even permutation of the 
Lorentz indices whereas £aj3-yS = — 1 for an odd permu- 
tation of the Lorentz indices and finally Sap-yS — oth- 
erwise. Second, the coefficients Bnj — B[Xi = A/ = 1) 
as well as T^nj — 'D{Xi = A/ = 1) contained various non 
contracted tensors. For example in there are thirty 
one non contracted tensors involving Sap-yS whereas in 
T^nj, there are sixty four. Particle physicists are very 
often dealing with these. As for us, we were confronted 
for the first time with coefficients like for example in Bnj 



and 



£(ai, a2, k, v) = £ai3-t5aia2k^v' 



(69) 



(70) 



Such coefficients are, at first sight, very complicated to 
evaluate. But, when following the conventions of A.G 
Grozin we found that the only non- vanishing contri- 
bution of oi corresponds to the Lorentz index a = \ while 
the only non- vanishing contribution of 02 corresponds to 



o 



DCS(up)+DCS(down) 
Unpolarized DCS 




-150000-100000 -50000 50000 100000 150000 

Photon Energy Transfer 



FIG. 3: Envelope of the relativistic DCS scaled in 10~^ a.u 
as a function of the energy transfer scaled in units of the laser 
photon energy lj for an electrical field strength oi £ = 1 a.u 
and a relativistic parameter 7 = 2. 



the Lorentz index (3 — 2. These helpful values are due 
to the choice of the 4-potential . Also, it allows to 
deduce that the only non trivial indices remaining for A:^ 
and are the pairs (7 = 0, J = 3) and (7 = 3, ^ = 0). 
So, the non contracted tensor e{ai,a2,k,pi) reduces to : 



e{ai,a2,k,pi) = |ap [£1203^°?? + £1230^^^°] 



|a| - 

c 



|pi| cos6'i - — 
c 



(71) 



where 9i is the angle between the initial electron mo- 
mentum and the Oz axis. The coefficient £(ai, 02, fc, v) is 
easier to evaluate since the Lorentz index corresponding 
to V = (1,0,0,0) is zero 6 = 0. This trick is often used 
in all QED textbooks JJ. So, the only choice left for k is 
7 = 3. Thus, we have : 



e{ai,a2,k,v) = ei23Qalalk^ 



(72) 



It is not necessary to give the whole set of the coeffi- 
cients involving the non contracted tensors appearing in 
Bn.f and 2?„j. It is sufficient to follow the rules concern- 
ing the tensor Eajs-fS and the geometry chosen for ai , a2 
and k, bearing in mind that = laje^* = |a|(0, 1, 0, 0), 
= jale!^ = |a|(0, 0,1,0) and finally fc'^ = f (1,0,0, 1). 
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The same holds for the coefficients Bf = B{\i = —A/ = 
1) and Vf — T>{Xi — —Xj = 1). The complete expres- 
sions of the coefficients A, B, C and V in both cases 
(n/, /) can explicitly be found in the Appcndix.lt is im- 
portant to remind that these coefficients occur only in 
presence of a laser field. If we put ~ (0,0), we eas- 
ily recover the results of section 2. Now, let us turn 
to the discussion of the results concerning dcr/dri/ 1 ^ ^ 
and dcr/dri/ 1 . We chose two regimes : a) the rela- 
tivistic regime corresponding to a relativistic parameter 
7 = 2 and an electric field strength £ — 1 a.u and b) 
the non relativistic regime corresponding to a relativis- 
tic parameter 7 = 1.0053 and an electric field strength 
£ = 0.05 a.u. There are two consistency checks that 
must be made. First, we have to show numerically that 
the sum of (DCS) / and {DCS)nj always gives the unpo- 
larized DCS and second, the non relativistic description 
for the unpolarized DCS must give the unpolarized DCS 
we have found in the formalism we have developed. 



A. The relativistic regime 



1,0 -, 




-0,2 -t-^ — I — ■ — I — ■ — I — ■ — I — ■ — I — ■ — I 

-180 -120 -60 60 120 180 



9, 



FIG. 4: Degree of polarization P for an electrical field 
strength of £ = 1. a.u and a relativistic parameter 7 = 2. 



This regime corresponds to a total energy of the incoming 
electron Ei — 0.510858 MeV which is the rest energy of 
the electron so, it is located below the threshold energy 
needed for the process of pair creation. As for the cutoffs. 



the negative part of the spectrum corresponds approxi- 
mately to —150000 photons (emitted by the laser). As a 
side-result, the various DCS (n, / or /) decrease abruptly 
because the arguments of the Bessel functions are close to 
their indices. Once again, there is an asymmetry between 
the emission part of the envelope and the absorption part 
of the same envelope for the unpolarized DCSs because 
of the presence of the denominator |q|'' in Ea. (p7|) . Now, 
the positive part of the spectrum corresponds approxi- 
mately to 175000 photons absorbed by the laser field. In 
Fig. (3), we show the envelope of da/(Klf as a function 
of the net number of photons exchanged. The cutoffs are 
n ~ —150000 photons for the negative part of the enve- 
lope and Umax — 170000 photons for the positive part 
of the envelope that corresponds to the absorptive part 
of the spectrum. The geometry chosen for Fig. (3) is 
Oi = 45°, 0, = 0°, 6f = 50° and 0/ = 90°. One many 
ask legitimately if such cutoffs are geometry-dependent. 
The answer is that indeed they are geometry dependent 
since for a different choice of the initial, final momen- 
tum angular coordinates Pi and p/ namely 9i = 30°, 
0j = 0°, Oi = 75° and 0/ = 90° the cutoffs obtained 
are rimin — —147000 photons and rimax — -1-154000 pho- 
tons. In Fig. (4), we show the degree of polarization P 
corresponding to this regime. 



DCS(up) 

DCS(down) 

5 r Unpolarized DCS 




-180 -120 -60 60 120 180 



FIG. 5: The various relativistic DCSs scaled in 10^^^ a.u as 
a function of the angle Of'm degrees for an electrical field 
strength oi £ = 1. a.u and a relativistic parameter 7 = 2. 
The corresponding number of photons exchanged is ±500. 
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DCS(up)+DCS(down) 
Unpolarized DCS 




a very important consistency check. Both Figs. (5,6) are 
accurate illustrations of this consistency check. 



B. The non relativistic regime. 

In this regime, the dressing of the angular coordinates of 
Pi as well as those of p/ is not important. The second and 
final consistency check is shown in Fig. (7). First, we give 
the non relativistic DCS obtained by using Schrodinger- 
Volkov wave functions : 



_ 4Z^|p/| 
dflf I Pi I 



1 



- p/ +nk|4 
(74) 

The two DCSs are very close, which was to be expected 
but with small deviations coming from the unpolarized 
relativistic DCS in the limit of small velocities and this is 
due to the fact that error propagation is more likely to oc- 
cur when one is dealing with very long expressions. This 
is indeed the case for the latter. We have summed over 
±500 photons to obtain the curves in Fig. (7). There is a 



FIG. 6: First consistency check related to the sum of DCS(up) 
and DCS (down) compared to the unpolarized DCS scaled in 
10~^^ a.u as a function of the angle Of in degrees for an elec- 
trical field strength of £^ = 1. a.u and a relativistic parameter 
7 = 2. The corresponding number of photons exchanged is 
±500. 



For this degree of polarization, we obtained a qualita- 
tive result similar to that of a linearly polarized field. We 
reached the same conclusion as in [9!| that the degree of 
polarization. 



[d<j/dnf{^) + da/dnf{i)] ' ' 

is weakly dependent on the number of photons ex- 
changed. This degree of polarization P varies as a func- 
tion of the angle df for the following geometry (Qi — 45°, 
(/)^ = 0, 0^ = 90° and < 6*/ < 180°). To begin with, 
we have made simulations concerning the various DCSs 
for a set of net number of photons exchanged. These 
sets (±100, ±200, ±300, ±400, ±500, ±1000) showed 
that the order of magnitude of the non flip DCS is close 
to the unpolarized DCS but the contribution of the flip 
DCS is not completely negligible even if it is small com- 
pared to both. The behavior of the three DCSs when the 
number of photons exchanged increases has an influence 
over the numerical values of the differential cross sections 
but as we are limited in our computational capabilities, it 
is not possible to achieve numerical convergence but the 
most important result is that the sum of DCS (up) and 
DCS (down) always gives the unpolarized DCS which is 



Non Relativistic DCS 

Relativistic DCS for small velocities 



3,0- 




0,0-1 — , — ^ — , — ^ — . — I — . — ^ — , — ^ — , — I 

-180 -120 -60 60 120 180 

9, 



FIG. 7: Second consistency check in the non relativistic 
regime related to the sum of DCS(up) and DCS(down) com- 
pared to the non relativistic DCS scaled in 10"* a.u as a func- 
tion of the angle in degrees for an electrical field strength 
oi E = 0.05 a.u and a relativistic parameter 7 = 1.0053. The 
corresponding number of photons exchanged is ±500. For 
large negative and positive angles, it becomes scaled in 10"'' 
a.u. 
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very good argument between the two curves and it is pos- 
sible by increasing the number of net photons exchanged 
to reach the same value as that of the unpolarized DCS 
in absence of the laser field. Many similar results were 
obtained because of a pseudo sum-rule that was shown 
by Bunkin and Fcdorov as well as KroU and Watson [l3] ■ 



3,5 -| 
3,0- 




Photon Energy Transfer 



FIG. 8: Envelope of the two unpolarized DCS in the non rel- 
ativistic regime for a number of net photons exchanged ±500 
scaled in units of lO"'^. The two curves are indistinguishable 
and are shifted from the elastic peak corresponding to 6'/ = 0. 
The cutoffs are visible. 

Even if the non relativistic description does not take 
into account the spin of the electron, the agreement be- 
tween the two DCSs shows in the one hand, that spin ef- 
fects are not important in this regime while on the other 
hand, given the more complete and sophisticated rela- 
tivistic description, the non relativistic limit of the rel- 
ativistic DCS always gives results very close to the non 



relativistic description. This does not mean that spin ef- 
fects are irrelevant but only that their contribution to the 
relativistic DCS (for small velocities) is small enough to 
be noticeable. As for the cutoffs, one can observe from 
Fig. (8) that for the negative part of the spectrum, it is 
located approximately at —450 photons while for the ab- 
sorptive part of the spectrum, it is located approximately 
at -1-470 photons. 

V. CONCLUSION. 



In this work, we have checked that results already known 
in the absence of a laser field are also valid in presence 
of a laser field, at least in first order of perturbation the- 
ory. The motivation of this study was to compare the 
order of magnitude for the various DCSs in the relativis- 
tic regime that we already found in [Q] particularly in the 
case of a laser field with linear polarization and we have 
found (see Fig. (6)) that indeed when scaled in 10*^~^^\ 
both DCSs are of similar order of magnitude and shape. 
Moreover, when focusing only on the formalism for po- 
larized electrons, the main difference between the linear 
and the circular polarization of the laser field, is that a 
thorough analysis of the contents of the formal expres- 
sion of DCS (up) and DCS (down) is more complicated in 
the case of the circular polarization than the linear po- 
larization of the laser field and this is due to the fact that 
these two DCSs contain non contracted symbols that we 
had to deal with for the first time. We succeeded in the 
simulations of the main two consistency checks namely 
that the sum of DCS(up) and DCS(down) always gives 
the unpolarized DCS regardless of the number of net pho- 
tons exchanged and also that in the non relativistic limit, 
both unpolarized DCSs (non relativistic and relativistic) 
give very close results. We also gave the envelopes of the 
two unpolarized DCSs in both regimes for the sake of il- 
lustration and also to retrieve the shape of Fig. (3) of (Tl| 
which was confirmed since the scaling is the same (10""^) 
and the asymmetric deviation from the elastic peak was 
also found as expected. The envelope for the relativistic 
regime was more difficult to obtain since the stability of 
the very robust code that evaluates the ordinary Bessel 
functions was put on trial because of the very high or 
very low arguments of the Bessel functions. Needless to 
say that working in simple precision, the results obtained 
are nonetheless sound and coherent. 



VI. APPENDIX 



A{Xi, Xf) = i^2.O*{k.pf)'^*{k.pi)'^*Xf*\i*\pf\'^*\pi\'^*c^*cos{0if)-2.O*{k.pf)'^*{k.p,)'^*\f*Xi*\pf^^ 
2.0*(/c.p/)^*(fc.pi)^*A/*Ai*|p/|*|pi|*c^'^ — 2.0*(fc.p/)^*(fc.pi)^*A/*Ai*|pip*c^*cos(^i/)*i?|-l-2.0^ 
Ai*c^*cos(0i/)*i?/*ii^i-l-2.O(/c.p/)^*(fc.pi)^*A/*Ai*c''*cos(6'i/)*i!J|*£^f *-|-2.0*(fc.p/)^* 
2.O*(fc.p/)^*(fc.pi)^*c-^^-)-2.O(fc.p/)^*(fc.pi)^*c**£'/*£'i + 2.O*(fc.p/)^*(fc.pi)*A/*Ai*(a^)*|pip*c'**cos(0^ 
2.O*(fc.p/)^*(fc.pi)*A/*Ai*(a^)*c^*cos(0i/)*i<^/*i?,j^*w--2.O*(fc.p/)^*(fc.pi)*(a^)*c^*£'i*w-l-2.O* 
Ai*(a^)*|p/P*c^*cos(0i/)*£'i*a; — 2.0*(fc.pj)*(fc.pi)^>i=A/*Ai*(a^)*c^*cos(6'i/)*i??*£'i*a; — ^ 
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c^*Ef*u — 2.0* {k.pf)*{k. Pi)* {k.Sf)*Xf* Xi*{a^)*\pf\*\pi\'^*c^*cos{9if)*iJ + 2.0*{k.pf)*{k. Pi)* {k.Sf)*Xf* Xi*{a^ 
\pf\*c^*cos{9if)*Ef*uj—2.0*{k.pf)*{k.pi)*{k.Sf)*Xf*Xi*{a'^)*\pi\*c^*uj—2.0*{k.pf)*{k.pi)*{k.Si)*Xf*Xi*{a'^)*\pf\'^* 
\pi\*c^*cos{9if)*uj — 2.0*{k.pf)*{k.pi)*{k.Si)*Xf*Xi*{a'^)*\pf\*c^*uj + 2.0*{k.pf)*{k.pi)*{k.Si)*Xf*Xi*{a'^)*\pi\*c'^* 
cos{9if)*Ej*uj — {k.pf)*{k.pi)*Xf*Xi*{a'^)'^*\pf\*\pi\*c^*i^'^ + (k.pf)*{k.pi)*Xf*Xi*{a'^)'^*cos{6if)*Ef*Ei*u^ 
{k.pf)*{k.pi)*Xf*Xi*{a^)*\pf\^*\pi\^*c'^*cos{9if)*(jj'^ — 2.0*{k.pf){k.pi)*Xf*Xi*{a^)*\pf\^*c^*cos{9if)*Ef*co'^ + 2.Q* 
{k.pf){k.pi)*Xf*Xi*{a^)*\pf\*\pi\*c^*Lj'^—2.0*{k.pf)*{k.pi)*Xf*Xi*{a'^)*\pi\'^*c^^ 

Xf*Xi*{a'^)*c'^*cos{9if)*Ef*Ei*ij^ + 2.0*{k.pf){k.pi)Xf*Xi*{a'^)*cos{9if)*Ej*Ef*iv'^ + {k.pf){k.pi)*{a'^)'^*c^*ijj'^ — 
2.0*{k.pf)*{k.pi)*{a'^)*\pf\*\Pi\*c^*cos{0if)*u'^ — 2.0{k.pf){k.pi)*{a'^)*c^*Lo^+2.0*{k.pf)*{k.pi)*{a'^)c'^*Ef*Ei*u^ 
{k.pf)*{k.Si)*Xf*Xi*{a'^) *\pf\*c'^*Ei*uj'^ — {k.pf){k.Si)*Xf*Xi*{a'^)'^*\Pi\*c'^*cos{9 

{a'^)'^*\pf\*(^*coa{9if)*Ei*iv'^ + {k.pi)*{k.sj)*Xf*Xi*{a'^)'^*\pi\*c'^*Ef*u}'^ + {k.sj)*{k.Si)Xf*Xi*{a'^)'^*\^ 
cos{9if)*Lo'^ + {k.Sf){k.Si)*Xf*Xi*{a'^)'^*c^*Lo'^ — {k.Sf){k.Si)*Xf*Xi*{a'^)'^*c'^*Ef*Ei*Lo'^^/[2.0*{k.pf)'^*{^ 

B{Xi,Xf) = ^{ai.pf)*{ai.pi)*{k.pf)'^*{k.Sf)*{k.Si)*Xf*Xi*c^ — 2.0*{ai.pf)*{ai.pi)*{k.pf)'^*{k.Si)*Xf*Xi*\pf\^ 

u-2.0*{ai:pf)*{ai:pi)*{k.pj)*{k.pi)*{k.Sf)*{k.Si)*Xf*Xi*c^+2M{ai.pf)*{ai.pi)*{k.pj)*{k.pi)*{k.sj)*Xf*X,*^^ 
c^*u) + 2.0*{ai.pf)*{ai.pi)*{k.pf)*{k.pi)*{k.Si)*Xf*Xi*\pf\*c^*uj — 2.0*{ai.pf)*{ai.pi)*{k.pf)*{k.pi)*Xf*Xi*\pf\* 
\pi\*c'^*ui^+2.0*{ai.pf)*{ai.pi)*{k.pf)*{k.pi)*Xf*Xi*c^*cos{9if)*Ef*Ei*ui^+2.0*{ai.pf)*{ai.pi)*{k.pf)*{k.pi)*c^* 
uj'^ + {ai.pf)*{ai/pi)*{k/pi)'^*{k.Sf)*{k.Si)*Xf*Xi*c^ — 2.0*{ai.pf)*{ai.pi)*{k.pi)'^*{k.Sf)*Xf*X 

{ai. Si)* {k. pf )^ *{k.pi)*{k.Sf)*Xf*Xi*c^ +2.0* {ai.pf)*{ai. Si)* {k.pf)'^*{k.pi)*Xf*Xi*\pf\*c^*LO+2.0*{ai.pf)*{ai. Si)* 
{k.pf)*{k.pi)'^*{k.Sf)*Xf*Xi*c^ — 2.0*{ai.pf)*{ai.Si)*{k.pf)*{k.pi)'^*Xf*Xi*\pf\*c^*w — 2.0* {ai.pf)*{ai.Si)^ 
{k.pi)*{k.sj)*Xj*Xi*c^*Ei*u}+2.0*{ai.pf)*{ai.Si)*{k.pf)*{k.pi)*Xf*Xi*\pf\*c'^*Ei*uj'^ — 2.0*{ai.pf)*{ai.Si)*{k.pf)* 
{k.pi)*Xf*Xi*\pi\*c^*cos{9if)*Ef*u!'^ — {ai.pf)*{ai.Si)*{k.pi)^*{k.Sf)*Xf*Xi*c^+2.0*{ai.pf)*{ai.Si)*{k.pi)'^*{k.Sf)* 
Xj*Xi*c^*Ei*uj—{ai.pi)*{ai.Sf)*{k.pf)^*{k.Si)*Xf*Xi*c^ + 2.0*{ai.pi)*{ai.Sf)*{k.pf)^*{k.pi)*{k.Si)*Xf*Xi*c^ — 
2.0*{ai.pi)*{ai.Sf)*{k.pf)'^*{k.pi)*Xf*Xi*\pi\*c^*uj + 2.0*{ai.pi)*{ai.Sf)*{k.pf)^*{k.Si)*Xf*Xi*c^*Ef*uj—{ai.pi)* 
{ai.Sf)*{k.pf)*{k.pi)'^*{k.Si)*Xf*Xi*c^+2.0*{ai.pi)*{ai.Sf)*{k.pf)*{k.pi)'^*Xf*Xi*\pi\*c^*u} — 2.0*{ai.pi)*{ai.Sf)* 
{k.pf)*{k.pi)*{k.Si)*Xf*Xi*c^*Ef*LLj~2.0*{ai.pi)*{ai.Sf)*{k.pf)*{k.pi)*Xf*Xi*\pf\*c'^*cos{9if)*Ei*uj^+2.0*{ai.pi)* 
{ai.Sf)*{k.pf)*{k.pi)*Xf*Xi*\pi\*c'^*Ef*u}^ + {ai.Sf)*{ai.Si)*{k.pf)^*{k.pi)*Xf*Xi*c^ — 2.0*{ai.Sf)*{ai.Si)*{k.pf)^* 
{k.pi)'^*Xf *Xi*c^ — 2.0* {ai.Sf)*{ai. Si)* {k.pf)'^*{k. Pi)* Xf*Xi*c^*Ef *u} + 2.0*{ai.s f)*{ai. Si)* {k.pf)'^*{k. Pi)* Xf* 
Xi*c^*Ei*Lu + (ai.Sf)*{ai.Si)*{k.pf)*{k.pi)^*Xf*Xi*c^ + 2.0*{ai.Sf)*{ai.Si)*{k.pf)*{k.pi)'^*Xf*Xi*c^*Ef*uj — 2.0* 
{ai.Sf)*{ai.Si)*{k.pf)*{k.pi)^*Xf*Xi*c^*Ei*Lu+2.0*{ai.Sf)*{ai.Si)*{k.pf)*{k.pi)*Xf*Xi*\pf\*\pi\*c^*cos{9if)*Lo^ + 
2.0*{ai.Sf)*{ai.Si)*{k.pf)*{k.pi)*Xf*Xi*c^*u'^ — 2.0*{ai.Sf)*{ai.Si)*{k.pf)*{k.pi)*Xf*Xi*c'^*Ef*Ei*u)'^ + {a2.pf)* 
{a2.pi)*{k.pf)^*{k.Sf)*{k.Si)*Xf*Xi*c^ — 2.0*{a2.pf)*{a.2-Pi)*{k.pf)'^*{k.Si)*Xf*Xi*\pf\*c^*LU—2.0*{a2.pf)*{a2.Pi)* 
(k.pf)* {k. Pi)* {k. s f)* {k. Si)* Xf*Xi*c^ + 2.0* {a2.pf)*{a2.pi)*{k. pf)*{k. Pi)* {k.s f)*Xf*Xi*\pi\*c^ *u} + 2.0* {a2.pf)* 
{a2-Pi)*ik.pf)*{k.pi)*{k.Si)*Xj*Xi*\pj\*c^*uj~2.0*{a2.pf)*ia2.pi)*{k.pf)*{k.pi)*Xf*Xi*\pf\*\pi\*c'^*uj'^+2.0*{a2.pf)* 
{a2.Pi)*{k. pf)*{k. Pi)* Xf*Xi*c^*cos{9if)*Ef*Ei*ijj'^ + 2.0* {a2.pf)*{a2.pi)*{k.pf)*{k.pi)*c^*u}^ + {a2.pf)*{a2.pi)* 
{k.pi)'^*{k.Sf)*{k.Si)*Xf*Xi*c^—2.0*{a2.pf)*{a2.pi)*{k.pi)'^*{k.Sf)*Xf*Xi*\pi\*c^*uj—{a2.pf)*{a2.Si)*{k.pf)'^*{k.pi)* 
{k.Sf)*Xf*Xi*c^+2.0*{a2.pf)*{a2.Si)*{k.pf)^*{k.pi)*Xf*Xi*\pf\*c^*iv+2.0*{a2.pf)*{a2.Si)*{k.pf)*{k.pi)^*{k.Sf)* 
Xf *Xi*c^ — 2.0* {a2.pf)*{a2.Si)*{k.pf)*{k.pi)'^*Xf*Xi*\pf\*c^*u} — 2.0*{a2.pf)*{a2.Si)*{k.pf)*{k. Pi)* {k.Sf)*Xf* Xi* 
c^*Ei*u)+2.0*{a2.pf)*{a2.Si)*{k.pf)*{k.pi)*Xf*Xi*\pf\*c'^*Ei*ui'^ — 2.0* {a2.pf)*{a2.Si)*{k.pf)*{k. Pi)* Xf*Xi*\pi\* 
c^*cos{9if)*Ef*u!'^—{a2.pf)*{a2-Si)*{k.pi)^*{k.Sf)*Xf*Xi*c^+2.0*{a2.pf)*{a2.Si)*{k.pi)'^*{k.Sf)*Xf*Xi*c^*Ei*uj— 
{a2.Pi)*{a2.Sf)*{k.pf)^*{k.Si)*Xf*Xi*c^+2.0*{a2.pi)*{a2.Sf)*{k.pf)'^*{k.pi)*{k.Si)*Xf*Xi*c^ — 2.0*{a2.pi)*{a2.Sf)* 
{k.pf)'^*{k.pi)*Xf*Xi*\pi\*c^*cj+2.0*{a2.pi)*{a2-Sf)*{k.pf)'^*{k.Si)*Xf*Xi*c^*Ef*tJ—{a2.pi)*{a2.Sf^ 
{k.Si)*Xf*Xi*c^+2.0*{a2.pi)*{a2.Sf)*{k.pf)*{k.pi)'^*Xf*Xi*\pi\*c^*LO — 2.0*{a2.pi)*{a2.Sf)*{k.pf)*{k.pi)*{k.Si)*Xf* 
Xi*c^*Ef*LO—2.0*{a2.pi)*{a2.Sf)*{k.pf)*{k.pi)*Xf*Xi*\pf\*c^*cos{9if)*Ei*LL!^+2.0*{a2.pi)*{a2.Sf)*{k.pf)*{k.pi)* 
Xf*Xi*\pi\*c^*Ef*uj'^ + {a2.Sf)*{a2.Si)*{k.pf)'^*{k.pi)*Xf*Xi*c^ — 2.0*{a2.s f)*{a2.Si)*{k.pf)^*{k.pi)'^*Xf*Xi*c^— 2.0* 
{a2.Sf)*{a2.Si)*{k.pf)'^*{k.pi)*Xf*Xi*c^*Ef*uj+2.0*{a2.Sf)*{a2.Si)*{k.pf)'^*{k.pi)*Xf*Xi*c^*Ei*uj+{a2.Sf)*{a2.Si)* 
{k.pf)*{k.pi)^*Xf*Xi*c^+2.0*{a2.Sf)*{a2.Si)*{k.pf)*{k.pi)'^*Xf*Xi*c^*Ef*ui—2.0*{a2.Sf)*{a2.Si)*{k.pf)*{k.pi)'^* 
Xf *Xi*c^* Ei*uj + 2.0*{a2.Sf)*{a2.Si)*{k.pf)*{k. Pi)* Xf*Xi*\pf\*\pi\*c^*cos{9if)*uj'^ + 2.0* {a2.Sf)*{a2. Si)* (k.pf)* 
{k.pi)*Xf *Xi*c^*uj'^ — 2.0* {a2.Sf)*{a2. Si)* {k.pf)*{k.pi)*Xf *Xi*c^*Ef*Ei*uj'^ + {k.pf)^*{k.pi)*Xf*Xi*{a'^)*\pf\*\pi\* 
c^ + {k.pf)^*{k.pi)*Xf*Xt*{a'^)*c^*cos{9if)*Ef*Ei + {k.pf)^*{k.pt)*{a'^)*c^ — {k.pf)^*{k.Si)*Xf*Xi*{a'^)*\pf\*c^*Ei — 
{k.pf)'^*{k.Si)*Xf *Xi*{a'^) *\pi\*c^ *cos{9i f)*Ef — 2.0* {k.pf)^*{k.pi)^*Xf *Xi*{a'^)*\pf\*\pi\*c^ — 2.0*{k.pf)^*{k.pi)-^* 
Xf*Xi*{a'^)*c'^*cos{9if)*Ef*Ei — 2.0*{k.pf)'^*{k.pi)'^*{a'^)*c^ — {k.pf)'^*{k.pi)*{k.Sf)*Xf*Xi*{a^)*\pf\*c^*cos{9if)* 
Ei — {k.pf)'^*{k.pi)*{k.Sf)*Xf*Xi*{a'^)*\pi\*c^*Ef+2.0*{k.pj)'^*{k.pi)*{k.Si)*Xf*Xt*{a'^)*\pf\*c^*Ei+2.0*{k.pf)'^* 
{k.pi)*{k.Si)*Xj*Xi*{a'^)*\pi\*c^*cos{9if)*Ef +2.0* {k.pf)'^*{k.pi)*Xf*Xi*{a'^)*\pf\'^*c^*cos{9if)*Ei*uj—2.0*{k.pf)'^* 
{k.pi)*Xf*Xi*{a^)*\pi\^*c'^*cos{9if)*Ef*co—2.0*{k.pf)'^*{k.pi)*Xf*Xi*{a'^)*c'^*cos{9if)*Ej*Ei*Lu+2.0*{k.pf)'^*{k.pi)* 
Xf*Xi*{a'^)*c^*cos{9if)*Ef*Ef*u—{k.pf)'^*{k.pi)*Xf*{a'^)*\pf 1*0^*0)— {k.pf)'^*{k.pi)*Xf*{a'^)*c^*cos{9f)*Ef*u)+ 
{k.pf)'^*{k.pi)*Xi*{a'^)*\pi\*c'^ *Lu—{k.pf)'^*{k.pi)*Xi*{a'^)*c^*cos{9i)*Ei*u!—2.0*{k.pf)'^*{k.pi)*{a'^)*c^*Ef*ui+2.0* 
{k.pf)-^*{k.pi)*{a^)*c^*Ei*co+{k.pf)^*{k.Sf)*{k.Si)*Xf*Xi*{a^)*\pf\*\pi\*c^*cos{9if) + {k.pf)^*{k.Sf)*{k.Si)*Xf*Xi* 
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(a^) *c^° + (fc.p/)^ * (fc. s/) * (fc. Si) * A/ * Aj *(a^)*c^*i?/*i?i — 2.0* (fc.p/)^*(fc.Si)* A/* Ai^ 
2.0*(fc.p/)^*(A;.Sj)*A/*Aj*(a^)*|p/|*c^*w+2.0*(fc.pj)^*(A:.Sj)*A/*Aj*(a^)*|pj|*c^*cos 
{a'^)*\pf\*c^*cos{9f)*u+{k.pf)'^*{k.Si)*Xi*{a'^)*\Pi\*c^*cos{9i)*u)—{k.pf)'^*{k.Si)*Xi*{a^ 

Xi*{a'^)*c''' *Ei*iv+{k.pf)*{k.pi)^*Xf*Xi*{a^)*\pf\*\pi\*c^ + {k.pf)*{k.pi)^*Xf*Xi*{a'^)*c'^*cos^^ 
{k.pi)^*{a'^)*c^+2.0*{k.pf)*{k.pi)^*{k.Sf)*Xf*Xi*{a'^)*\pf\*c^*cos{6if)*Ei+2.0*{k.pf)*{k.pi)'^*{k.s^^^ 
\pi\*c^*Ef — {k.pf)*{k.pi)^*{k.Si)*Xf*Xi*{a'^)*\pf\*c^*Ei — {k.pf)*{k.pi)'^*{k.Si)*Xf*Xi*{a'^)*\pi\*c^*cos{^ 
2.0*{k.pf)*{k.pi)'^*Xf *Xi*(a'^)*\pf \'^*c**cos{6if) *Ei*oj +2.0* {k.pf)*{k. pi)'^*Xf*Xi* {a'^)*\pi\'^*c**m 
2.0*{k.pf)*{k.piY*Xf*Xi*{a'^)*c^*cos{9if)*E'j*Ei*u — 2.Q*{k.pf)*{k.piY*Xf*Xi*{a^)*c^*cos{^ 
{k.pf)*{k.piY*Xf*{a'^)*\pf\*c'*ui — 2.{]*{k.pf)*{k.piY*Xi*{a'^)*\pi\*c'*ui+2.{]* 
{k.pf)*{k.piY*{a^)*c^*Ei*u) — 2X]*{k.pj)*{k.pi)*{k.Sf)*{k.Si)*Xf*Xi*{a'^)*\pj\*\pi\*c^ 
{k.Sf)*(k.Si)*Xf*Xi*(a^)*c^'-' — 2.0*{k.pf)*{k.pi)*{k.Sf)*(k.Si)*Xf*Xi*(a^)*c^'*Ef*Ei+2.0*{k.pf)*{^ 
Xi*{a^)*\pf\*\pi\'^*c^''*cos{9if)*iv — 2.Q*(k.pf)*{k.pi)*{k.Sf)*Xf*Xi*(a'^)*\pf\*c**cos{9if)*Ef*u^ 
(k.Sf)*Xf*Xi*{a'^)*\pi\*c^*w + 2.0*{k.pf)*{k.pi)*{k.Si)*Xf*Xi*{a'^)*\pf\'^*\pi\*c^*cos{9if)*i^ 
{k.Si)*Xf*Xi*{a'^)*\pf\*c^*uj—2.0*{k.pf)*{k.pi)*{k.Si)*Xf*Xi*{a'^)*\pi\*c'^*cos{9if)*Ej*ij+2.^^ 
Xi*{a^)*c^ *Ef*Lu+2.0*{k.pf)*{k.pi)*(k.Si)*Xi*{a'^)*c^*Ei*ij—2.0*{k.pf)*{k.pi)*Xf*Xi*(a^ 
w^ + 2.0*(/c.p/)*(A:.pi)*A/*Ai*(a^)*|p/P*c^*cos(%)*i??*a;^ — 2.0*(fc.p/)*(fc.p,;)*A/*Ai^ 
{k.pf)*{k.pi)*Xf*Xi*{a'^)*\pi\'^*c^*cos{9if)*Ej*Lo'^+2.0*{k.pf)*{k.pi)*Xf*Xi*{a'^)*c'^*cos{9if)*Ef*Ei*u^ 

(A;.p/)*(A:.pi)*A/*Ai*(a^)*cos(6'i/)*i?|*i?f *cj^-2.0*(fc.p/)*(fc.pi)*A/*(a^)*|p/P*c^*cos(% 
Xf*{a'^)*\pf\*\pi\*c^*coa{9i)*iv'^ — 2.0*{k.pf)*(k.pi)*Xf*{a'^)*\pf\*c^*Ei*iv'^+2.0*{k.pf)*{k.pi)*Xf*{^^ 
Ej*iv^ — 2.0*{k.pf)*{k.pt)*Xi*{a'^)*\pi\*c^*Ef*Lv'^ + 2.0*{k.pf)*{k.pi)*Xt*{a'^)*c'^*cos{9t)*Ef*Et^ 
{k.pi)*{a'^)*\pf\*\pi\*c^*cos{9if)*u'^+2.0*(k.pf)*{k.pi)*(a'^)*c^*u)'^ — 2.0*{k.pf)*{k.pi^ 

{k.Sf)*Xf*Xi*{a'^)*\pf\*c^*cos{9if)*Ei — {k.pi)^*{k.Sf)*Xf*Xi*{a'^)*\pi\*c^*Ef — {k.pi)^*Xf*{a'^)*\pf\*c'^*iv+{k.^^ 
Xf*{a^)*c^*coa{9f)*Ef*u}+{k.pi)^*Xi*{a^)*\pi\*c^*oj+{k.pi)^*Xi*{a^)*c^*cos{9i)*Ei*iv+{k.pi)'^*{^ 
Xi*{a'^)*\pj]*\pi\*c^*cos{9if) + {k.pi)^*{k.Sf)*{k.Si)*Xf*Xi*{a'^)*c^" + {k.pi)^*{k.Sf)*{k.Si)*Xf*Xi*{a'^)* 
2.0* {k. pi)"^ * {k.s f)*Xf*Xi*{a'^)*\p f\*\pi\'^*c^ *cos{9if)*(jj+2. 0*{k.pi}'^*(k.s f)*Xf*Xi*(a'^)* 
2.0*{k.pi)'^*{k.Sf)*Xf*Xi*{a^)*\pi\*c^*oj+{k.pi)'^*{k.Si)*Xi*{a^)*\pf\*c^*cos{9f)*u — (k.pi)'^*(k.Si)*Xi*{^^ 
cos{9i)*iv—{k.pi)'^*{k.Si)*Xi*(a'^)*c'^*Ef*uj—{k.pi)'^*{k.Si)*Xi*(a'^)*c''^*Ei*ij+2.0*{k.pi)'^*Xf*{^^ 
2.0*{k.pi)^*Xf*{a'^)*c^*cos{9f)*Ef*Ei*Lj'^ — 2.0*{k.pi)^*Xi*{a'^)*\pf\*\pi\*c^*cos{9j)*Lu^+2.0*{k. 

c^*cos(6lj)*w2_^2.0*(fc.pj)2*Aj*(a2)*|pi|*c5*S/*tj2_2.0*(fc.p,)2*Ai*(a2)*c''*cos(^ 



C(Ai, Xf) = < cos(2(/)o)*(ai.p/)*(ai.pi)*(A:.p/)^*(fc.s/)*(A:.Si)*A/*Ai*c'^ — 2.0*cos(2(/)o)*(ai.p/)*(ai.pi)*(fc.p/)^* 



{k.Si)*Xf*Xi*\pf\*c^*ui—2.0*cos{2(j)o)*{ai.pf)*{ai.pi)*{k.pf)*{k.pi)*{k.Sf)*{k.Si)*Xf*Xi*c^+2.0*cos{2(j)o)*{ai.pf^ 
{a\.pi)*{k.pf)*{k.pi)*{k.Sf)*Xf*Xi*\pi\*d^*tJ+2.0*cos{2(l)o)*{ai.pf)*{ai.pi)*{k.pf)*{k.pi) 

2.0*cos{2(j>o)*iai.pf)*{ai.pi)*(k.pf)*{k.pi)*Xf*Xi*\pf\*\pi\*c^*uj'^+2.0*cos{24>o)*iai.pf)*{ai.pi)*{k.pf)*{k.p 
Xi*cos{9if)*Ef*Ei*oj'^+2.0*cos{2(j)o)*(n.i.pf)*(ai.pi)*{k.pf)*(k.pi)*c'^*(jj'^+cos{2(f)o)*{ai.pf)* /)* 
{k .Si) * X f * Xi*c^ — 2.0*cos{2(j)o)* {ai .p f) * {(1 1 .pi) * {k .pi)'^ * {k .3 f) * X f * Xi*\pi\*c^ — cos{2(po)* {ai -P f) * 
{k.pi)*{k.Sf)*Xf*Xi*c^+2.0*cos{2(j)o)*{ai.pf)*{ai.Si)*{k.pf)^*{k.pi)*Xf*Xi*\pf\*c^*ui+2.0*cos{2(j)o)*{ai.pf)*{ai.Si 
{k.pf)*{k.pi)'^*{k.Sf)*Xf *Xi*c^ — 2.0* cos{2(l)o)*{ai.pf)*{ai.Si)*{k.pf)*{k.pi)^*Xf*Xi*\pf\*c'^*LO — 2.0*cos{2(j)o)* 
{ai.pf)*{ai.3i)*{k.pj)*{k.pi)*{k.Sf)*Xf*Xi*c'^*Ei*LLj+2.0*cos{2(l)o)*{ai.pf)*{ai.Si)*{k.pf)*{k.pi)*Xf*Xi*\pf\*c^*E^^ 
uj"^ — 2.0*cos{2(j)o)*{ai.pf)*{ai.Si)*{k.pf)*{k.pi)*Xf*Xi*\pi\*c^*cos{9if)*Ef*u'^~cos{2(j)o)*{ai.pf)*{ai.Si)*{k.pi^^ 
{k.Sf)*Xf*Xi*c^+2.0*cos{2(l>o)*{ai.pf)*{ai.Si)*{k.pi)^*{k..Sf)*Xf*Xi*c'^*Ei*u — cos{2(j)())*{ai.pi)*{ai.s 
{k.Si)*Xf*Xi*c^+2.0*cos{2(j)o)*(ai.pi)*(ai.Sf)*(k.pf)'^*(k.pi)*(k.Si)*Xf*Xi*c^—2.0*cos{2(j)o)*{ai.pi)*{ai.Sf)*(k.pf^^^ 
{k.pi)*Xf*Xi*\pi\*c^*ui + 2.0*cos{2(po)*{ai.pi)*{ai.Sf)*{k.pf)^*{k.Si)*Xf*Xi*c^*Ef*Lv — cos{2(j)o)*{ai.pi)*{ai.Sf)* 
{k.pf)*{k.pi)^*{k.Si)*Xf*Xi*c^+2.0*cos{2(po)*{ai.pi)*{ai.Sf)*{k.pf)*{k.pi)'^*Xf*Xi*\pi\*c^*tJ—2.0*c^^ 
{ai.Sf)*{k.pf)*(k.pi)*{k.Si)*Xf*Xi*c^*Ef*uj—2.0*cos{2(j)o)*{ai.pi)*{ai.Sf)*{k.pf)*{k.pi)*Xf*Xi*\pf\*c'^*cos{9if)* 
Ei*iv'^+2.0*cos{2(po)*{ai.pi)*{ai.Sf)*{k.pi-)*{k.pi)*Xf*Xi*\pi\*c^*Ef*uj'^+cos{2(j>o)*{ai.Sf)*{ai.Si)*{k.pf)^*{k.pi)* 
Xf*Xi*c^ — 2.0*cos{2(j)o)*{ai.Sf)*{ai..Si)*{k.pf)^*{k.pi)^*Xf*Xi*c^ — 2.0*cos{2(j)o)*{ai.Sf)*{ai.Si)*{k.pf)^*{k.pi)* 
Xf*Xi*c^*Ef*uj + 2.0*cos{24>o)*{ai.Sf)*{ai.Si)*{k.pf)^*{k.pi)*Xf*Xi*c^*Ei*uj+cos{2(po)*{ai.Sf)*{ai.Si)*{k.pf)* 
{k.pi)^*Xf*Xi*c^+2.0*cos{2.0(l)o)*{ai.Sf)*{ai.Si)*{k.pf)*{k.pi)^*Xf*Xi*c'^*Ef*u) — 2.0*cos{2.0(j)())*{ai.Sf)*{ai.Si)* 
{k.pf)*{k.pi)^*Xf*Xi*c'^*Ei*u+2.0*cos{2.0*(j)o)*{ai.Sf)*{ai.Si)*{k.pf)*{k.pi)*Xf*Xi*\pf\*\pi\*c'^*cos{9if)*iv'^+2.0* 
cos{2(l)o)*{ai.Sf)*{ai.Si)*{k.pf)*{k.pi)*Xf*Xi*c^*uj'^ — 2.0*cos{2(j)o)*{ai.Sf)*{ai.Si)*{k.pf)*{k.pi)*Xf*Xi*c'^*Ef*Ei* 
u^ — cos{2(j)o)*{a2.pf)*{a2.pi)*{k.pf)'^*{k.Sf)*{k.Si)*Xf*Xi*c^+2.0*cos{2(j)o)*{a2.pf)*{a2.pi)*{k.pf)^*{k.Si)*Xf*Xi* 
\pf\*c**iv+2.0*cos{2(j)o)*(a2-Pf)*{<i2-Pi)*{k.pf)*{k.pi)*{k.Sf)*{k.Si)*Xf*Xi*c^~2.0*cos{2(po)*{a2.pf)*(a2.pi) 
{k.pi)*{k.Sf)*Xf*Xi*\pi\*c^*w—2.0*cos{2(l)o)*{a2-Pf)*{a2.pi)*{k.pf)*{k.pi)*{k.Si)*Xf*Xi*\pf\*c^*Lv+2.0*cos{2(l)o)* 
{a2-Pf)*{a2.pi)*{k.pf)*{k.pi)*Xf*Xi*\pf\*\pi\*c?*u)^—2.0*cos{2(l)o)*{a2-Pf)*{a2 .pi)*{k.pf)*{k.pi)*Xf*Xi*cos{9if)*Ef* 
Ei*Lu'^ — 2.0*cos{2(t>o)*{a2.pf)*{a2.pi)*{k.pf)*{k.pi)*c^*iv'^—coa{2(po)*{a2.pf)*(a2.pi)*{k.pi)'^*{k.Sf)*{k.Si)*Xf*Xi* 
c^+2.0*cos{2(j)o)*{a2-Pf)*{a2.pi)*{k.pi)^*{k.Sf)*Xf*Xi*\pi\*c'^*u!+cos{2(j)o)*{a2-Pf)*{a2.Si)*{k.pf)'^*{k.pi)*{k.Sf)* 
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A/*Ai*c^ — 2.0* cos(2</)o)*(a2. p/)*(a2. Si)* (A;.p/)^*(/c.pi)* A/*Ai*|p/| *c^*u; — 2.0* cos(2(/)o)*(a2.p/)*( 

(fc.Pi)^ * (fc. s j^) * A/ * Ai * + 2. * cos(2(/)o) * (a2 . j3/)*(o2. Si) *(fc.pj^)* (A: .pi)^* A/ *Ai*|p/ 1 *c^*aj + 2. 0*cos(2(/)o)*(a2.p/)* 

(a2-Si)*(fc.p/)*(fc.pi)*(fc.s /)*A f*Xi*c^*Ei*uj~2.0*cos{2<j>o)*{a2-Pf)*{a2-Si)*{k.pf)*{k.pi)*X f*Xi*\p j\*c^*Ei*uj'^+2.0* 

cos{2(po)*{a2-Pf)*{a2-Si)*{k.pf)*{k.pi)*Xf*Xi*\pi\*c'^*cos{0if)*Ef*w^+cos{2(j)o)*{a2.pf)*{a2-Si^ 

Ai*c^ — 2.0* cos(20o)*(a2-P/) *(a2 -Si)* (fc.pi)^*(A:.s/)* A/*Ai* c^*i5i*w +cos(2(^o)*(a2 -Pi)* (a2-s/)*(fc.p/)^* (fc. Si)* A/* 

Xi*c^ — 2.0*cos{2(/)o)*{a2.pi)*{a2-Sf)*{k.pf)'^*{k.pi)*{k.Si)*Xf*Xi*c^+2.0*cos{2(p())*{a2 

A/*Ai*|pi|*c^*w — 2.O*cos(2(/>o)*(a2.pi)*(a2-s/)*(fc-P/)^*(fc-Si)*A/*Ai*c^*i?/*a;+cos(20o)*(a2-Pi)*(a2-S/)*(fc.p/)* 
{k.pi)'^*{k.Si)*Xf*Xi*c^—2.0*cos{2(j)o)*{a2.pi)*{a2-Sf)*{k.pf)*{k.pi)^*Xf*Xi*\pi\*c'^*w 

{k.pf)*{k.pi)*{k.Si)*Xf*Xi*c'^*Ef*ij + 2.0*cos{2(f>Q)*{a2.pi)*{a2-Sf)*(k.pf)*{k.pi)*Xf*Xi*\pf\*c'^*cos{dif^ 
2.0*cos(2(/)o)*(a2 •Pi)*(a2.s/)*(A:.p/)*(fc.pj)*A/*Ai*|pi|*c^*-E'/*a;^ — cos(2(/>o)*(a2.s/)*(o2.Si)*(A;.p/)^*(fc.pi)*A/*Ai* 
c^+2.0*cos(2^o)*(oi2-S/)*(a2-Si)*(fc-i'/)^*(fc-ft)^*A/*Aj*c^+2.0*cos(2(^o)*(«2-S/)*(a2.Si)*(fc 

Ef — 2.0*cos{2(f>o)* {a2.s f)*{a2-Si)*{k.pf)'^*{k.pi)*Xf*Xi*c^*Ei*ij — cos{2(j)Q)* {a2-s f)*{a2 -Si)* (k.^ 

Ai*c^ — 2.O*cos(20o)*(a2-s/)*(a2.Si)*(fc.p/)*(fc.pi)-^*A/*Ai*c'**i?/*u; + 2.O*cos(20o)*(a2-S/)*(a2-Si)*(fc-P/)*(fc-Pi)^* 
A/*Ai*c''*ii^i*w— 2.O*cos(20o)*(a2-s/)*(o2.Si)*(A;.p/)*(A;.pi)*A/*Ai*|p/|*|pi|*c'**cos((?i/)*a;^ — 2.0*cos(2(/)o)*(a2.s/)* 
(a2-Si)*(fc.p/)*(fc.pi)*A/*Ai*c^*w^H-2.O*cos(20o)*(a2-s/)*(a2-Si)*(fc.p/)*(fc.pi)*A/*Ai*c^*i^/*£'i*w^+sin(2(/)o)* 
(ai. p/) * (a2 . Pi) * (fc.p/)^ * (A). s/) * (fc. Si) * A/ * Ai *c^ — 2.0* sin(2(^o)*(«i-P/)*(a2-Pi)*(fc-P/)^*(fc-Si)* A/* Ai*|p/ 1 *c^*w — 
2. * sin(2</)o) * (ai. p/) * (a2 . Pi) *(fc.p/)* (A: .pi)* (A;. s/)*(fc. Si)* A/ *Ai*c^ + 2. O*sin(20o)*(ai-P/)*(a2 -Pi) *(fc-?'/)*(fc-Pi)* 
(A:. s/) * A/ * Ai * |pi I *c'' *a>'+2. 0*sin(2(/)o) * (ai . p/)* (a2 . Pi) *(A;.p/)* (A: .pi)* (A;. Si )*A|*Ai*|p/ 1 *c^*a;— 2. 0*sin(2(/>o)*(ai-P/)* 
(a2.pi)*(A:.p/)*(A\pi)*Ay-*Ai*|pj|*|pi|*c^*w^ + 2.O*sin(20o)*(oi-P/)*(Q^2-Pi)*(fc-P/)*(^-Pi)*A/*Ai*cos(6'ij)*i?y*£Ji* 
+ 2. * sin(2(/)o) * (ai. p/ )* (a2 . Pi) * (A;.p/) * (A: .pi)*c^*cj^+sin(2(/)o)*(ai.p/)*(a2. Pi) *(A;.pi)^* (A:. s/)* (A:. Si)* A/ *Ai*c^ — 
2. 0* sin(20o) *(«! -P/) * ("2 -Pi) * ('c-Pi)^ * (fc-S/) * A/ * Ai*|pi I *c"'*w — sin(2^o)*(«i-P/)*(«2 -Si)* (A;.p/)^*(A;.pi)*(fc.s/)*A^ 
Ai*c^+2.0*sin(2(^o)*(«i-P/)*(a2-Sj)*(A;.p/)^*(fc.Pi)*A/*Aj*|p/|*c^*a;+2.0*sin(2(/)o)*(ai.p/)*(a2-Si)*(A;.p/)*(A;.pi 
(A:.s/)*A/*Ai*c^— 2. 0*sin(2(^o)*(ai •P/)*(a2 -Si)* (A: .p/)*(A;.pi)^* A/*Ai*|p/ |*c^*cj—2. 0*sin(2(/)o)*(ai.p/)*(a2. Si)* (A;.p/)* 

(A: .pi )* (A;. s/) * A/ * Ai * c'' * i?i * + 2. * sin(2(/)o )* (oi . p/) * (a2. Si )* (A: .py- )* (A; .pi)*A/* A; *|p| I *c^*_Ei*cj^ — 2.0* siii(20o)* 
(ai.p/)*(a2.Si)*(A:.p/)*(A;.pi)*A/*Ai*|pi|*c^*cos(0i/)*ii^/*w^ — sin(2(/)o)*(ai.p/)*(a2-Si)*(A;.pi)''^*(A:.s/)*A/*Ai*c^ + 
2. * sin(2(/)o) * (ai. p/) * (a2. Si) * (A:.pi)^ * (A:. s/) * A/ * Ai *c'^*i?i*a;+sin(2(/>o)*(ai -Pi) *(a2-P/)*(A:.p/)^* (A:. s/)* (A;. Si)* A/* 
Ai*c^ — 2.0 *sin(2(/)o)*(ai . Pi)* (a2.p/)*(A;.p/)^ * (A:. Si)* A/* Ai*|p/ 1* c^*w — 2.0 *sin(2(^o)*(ffli -Pi)* (a2-P/)*(A:.p/) *(A:.pi)* 
(A;. s/) * (A;. Si) * A/ * Ai *c^ +2. O*sin(20o) * (ai -Pi) * (^2 •P/)*(fc-P/)*(A:.pi)* (A;. s/)* A/ * Aj *|pi I *c'' *Li;+2.0*sin(2^o)*(ai -Pi)* 
(a2.p/)*(A;.p/)*(A;.p,;)*(fc.Si)*A/*Ai*|p/|*c^*a;— 2. 0*sin(2(/)o)*(ai . Pi)* (a2.p/)*(A;.p/)* (A: .pi)* A/* Ai*|p/ 1* |pi I* 0^*6^^ + 

2.0*sin(2(/)o)*(ai.pi)*(a2 J''/)*(fc-P/ )* (A: .pi)* A/* Ai*cos(0i/) *i?/*£'i* +2. O*sin(20o)*(ai -Pi)* (a2-P/)*(fc-P/) *(fc-Pi)* 
*w^+sin(2(/)o) * (ai. Pi) * (a2 . p/ )* (A: .pi)^ * (A:. s/ )* (A;. Si )* A/ * Ai *c^— 2. 0*sin(2(/)o)*(ai . Pi)* (a2.p/)* (A: .pi)^* (A:. s/)* A/* Ai* 
|Pi|*c^*w— sin(2^o)*(ai. Pi)* (a2-s/)*(A;.p/)^* (A;. Si)* A/* Ai*c^+2.0*sin(2^o)*(ai -Pi)* (02 •s/)*(A;.p/)^*(fc.pi)* (A;. Si)* A/* 
Ai*c^ — 2. 0*sin(2(^o)*(ai •pO*('^2-s/)* (A;.p/)^*(A;.pi)* A/* Aj* |pi I* c^* w+2.0*sin(2(/>o)*(ai. Pi)* (a2.s/)*(A;.p/)^* (A;. Si)* 
A/*Ai*c'**i?/*cj — sin(2^o)*(cn -Pi)* (a2 -s/)* (A: .p/)*(A;.pi)^* (A:. Si)* A/* Ai*c^ + 2.0* sin(2(/)o)*(ai. Pi)* (a2.s/)* (A: .p/)* 

(A:.pi)^*A/*Ai*|pi|*c''*a;— 2. O*sin(20o)*(ai -Pi)* (02 -s/)* (A: .p/)*(A;.pi)* (A;. Si)* A/* Ai*c^*i?/ *cj — 2. O*siii(20o)*(ai -Pi)* 
(a2. s/ )* (A: .p/ )* (A;.pi )* A/ * Ai * |p/ 1 * * cos(f?i/ )* i?i *<x'^ +2. * sin(20o)*(ai -Pi )*(a2-S/)* (A: .p/)*(A;.pi)* A/ *Ai*|pi I *c^* 
i?/*a;^ — sin(2(/)o)*(ai. s/)*(a2. Pi)* (A:.p/)"^* (A:. Si)* A/* Ai*c^ +2.0* sin(2(/)o)*(ai. s/)*(a2. Pi)* (A: .p/)^*(A:.pi)* (A:. Si)* A/* 
Ai*c^ — 2.0* sin(20o)*(ai •s/)*(a2. Pi)* (A:.p/)^*(A:.pi)* A/* Ai*|pi I* c^*w +2. O*sin(20o)*(«i -s/)* (02 . Pi)* (A:.p/)^* (A;. Si)* 
A/*Ai*c'**i?y*a; — sin(2(/)o)*(ai. s/)* (02. Pi)* (A: .p/)* (A; .pi)^* (A;. Si)* Ay* Ai*c^ + 2.0* sin(2(/)o)*(ai. s/)*(a2. Pi)* (A; .py)* 
(A:.pi)'^*A/*Ai*|pi|*c**cj — 2. O*sin(20o)*(ai •s/)*(a2.Pi) *(A;.p/)*(A;.pi) * (A;. Si)* Ay* Ai*c^*ii^/*a; — 2. 0*sin(2(/>o)*(ai.s/)* 
(a2.pi)*(A:.py)*(A:.pi)*Ay*Ai*|py|*c^*cos(6'iy)*£'i*a;^ + 2.0*sin(2(/)o)*(ai.sy)*(a2.pi)*(A:.py)*(fc.pi)*Ay*Ai*|pi|*c^* 
i?y *cj^ + sin(2(/)o)*(ai. sy)*(a2. Si)* (A;.py)'^*(A;.pi) * Ay *Ai*c^ — 2.0* sin(20o)*(ai •S/)*(a2. Si)* (A: .py)^* (A: .pi)^* A/* Ai* 
c^ — 2.0* sin(2^o)*(ai -s/)* (a2 -Si)* (A;.p/)^*(A;.pi)* A/* Ai*c''*ii^/*a; +2.0* sin(2^o)*(«i-s/) *(a2 -Si)* (A;.p/)^*(A;.pi)* A/* 
Ai*c**£'i*w+sin(2^o)*(ai-s/)*(a2-Si)*(A;.p/)*(A;.pi)^*A/*Ai*c® + 2.0*sin(2(^o)*((ii-s/)*(a2-Sj)*(A;.p/)*(A;.pi)^*A/* 
Ai*c'''*i?/*w — 2.0* sin(2(^o)*(ai -sy)* (a2 -Si)* (A: .py)* (A: .pi)^* Ay* Ai*c'**i?i *cl' + 2.0* sin(20o) *(ai-sy)* (a2 -Si)* (A: .py)* 
(A: .pi )* Ay * Ai * |py I * |pi I * c'' * cos(6'iy )* cj^ + 2. * sin(2(/)o) * (fli. sy )* (a2. Si)* (A: .py)* (A: .pi)* A/ *Ai*c^*a->^ — 2.0* sin(2(/)o)* 
(oi.s y )* (0,2. Si) * (A;.py)* (A;.pi )* Ay * Ai*c^* Sy * i?i* w^ — sin(2(/)o)*(ai. Si)* (a2.py)* (A- .py)^*(A;.pi)* (A:. s/)* Ay* Ai*c^ +2.0* 
sin(20o) *(«! -Si) * (02 . p/) * (A:.p/)^ * (A:.pi) * A/ * Ai *|p/ 1 * *a; +2. 0*sin(2(/i>o)*(ai. Si)* (02 .p/)*(A;.p/)*(A;.pi)^* (A;. s/)* A/* 
Ai*c^ — 2.O*sin(20o)*(ai-Si)*(«2-P/)*(A:.py)*(A;.pi)^*A/*Ai*|py |*c^*w — 2.O*sin(20o)*(«i-Si)*(a2-P/)*(A:.py)*(A;.pi)* 
(A;.s/)*Ay *Ai* c''*i?i* CJ + 2.0* sin(2(/)o)*(oi. Si)* (02. p/)* (A: .p/)* (A; .pi)* A/*Ai*|p/ 1* c^*_Ei*w^ — 2.0* sin(2(?!)o)* (fli. Si)* 
(a2.p/)*(A;.py)*(A:.pi)*A/*Ai*|pi|*c^*cos(6'iy)*i5y*a;^ — sin(2(/)o)*(ai. Si)* (02. p/)*(A;.pi)'^* (A;. sy)* A/* Ai*c^ + 2.0* 
sin(20o)*(ai-Si)*(a2-p/)*(fc-pi)^*(fc-sy)*Ay * Aj* c^*i^i *a;+sin(2(/>o)* (fli -Si)* (02. sy)* (A: .py)^* (A: .pi)* A/* Ai*c® — 2.0* 
sin(20o)*(ai-Si)*(a2-sy)*(A:.py)^*(A;.pi)-^*Ay*Ai*c^ — 2.0* sin(2(/)o)*(ai. Si)* (02. sy)*(A;.py)^ *(A;.pi)* Ay* Ai*c'**i?y *u; + 
2. 0*sin(2(/)o)*(ai -Si) *(a2-S/)*(A;.p/)^*(A;.pi)* A/ *Ai*c''*i?i*a;+sin(2^o)*(ai -Si)* (a2.s/)*(A;.p/)*(A;.pi)^* A/ *Ai*c^+2.0* 
sin(2(^o)*(oi-Si)*(a2-s/)*(A;.p/)*(A;.pi)^*A/*Ai*c**i^/*a;— 2.0*sin(2(^o)*(fli-Si)*(a2-s/)*(A;.p/)*(A;.pi)^*A/*Ai*c^*£'i* 
u;+2.0*sin(2(/)o)*(ai. Si )*(a2.s y)*(A:.py)*(A;.pi)*Ay*Ai*|py|*|pi|*c'**cos(f?iy)*cj^+2.0*sin(2(/)o)*(ai.Si)*(a2.s/)*(A;.py)* 

(A;.pi)*A/*Ai*c^*w^ — 2.0* sin(20o)*(«i -Si)* (fl2-s/)*(A;.p/)* (A;.pi)* A/* Ai*c^*i?/ *i?i*w^ W[2.0*(A;.p/)^*(A;.pi)^*c^] 
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I>(Ai, A/) = I cos{<po)*{ai.pf)*{k.pf)'^*{k.pi)'^*Xf*Xi*\pf\*\pi\*c^ + cos{<po)*{ai.pf)*{k.pf)'^*{k.pi)'^*Xf*\i*c'^* 

cos{9if)*Ef*Ei+cos{(po)*{ai.pf)*{k.pf)'^*{k.pi)'^*c^ — coa{(po)*{ai.pf)*{k.pf)'^*{k.pi)*{k.Si)*Xf*\i*\pf\*c^*Ei~ 

cos{(po)*{ai.pf)*{k.pf)'^*{k.pi)*{k.Si)*Xf*Xi*\pi\*c^*cos{9if)*Ef-cos{(po)*{ai.pf)*{k.pf)*{k.pi)^*Xf^ 

c^ — cos{(f>())*{ai.pf)*{k.pf)*{k.pi)^*Xf*Xi*c'^*cos{0if)*Ef*Ei — cos{(f)Q)*{ai.pf)*{k.pf)*{k.pi)^*c^+cos{(f)Q)* 

{k.pf)*{k.pi)'^*{k.Si)*Xf * Xi*\pf\*c^* Ei + cos{(j)o)* {ai.pf)*{k.pf)* {k.pi)'^*{k. Si)* Xf*Xi*\pi\*c^*coa{9if)*Ef — 2.0* 

cos{(j)o)*{ai.pf)*{k.pf)*{k.pi)'^*Xf*Xi*\pi\'^*c'^*cos{9if)*Ef*w + 2.0*cos{(j)o)*{ai.pf)*{k.pf)*{k.^^^ 

coa{9if)*Ef*Ef*Lj+2.0*cos{(j)o)*{ai.pf)*{k.pf)*{k.pi)'^*c^*Ei*ij—cos{(f)o)*{ai.pf)*{k.pf)*{k.pi)*{k.Sf) 

\Pi\*c'^*ij+cos{(j)o)*{ai.pf)*{k.pf)*{k.pi)*Xf*Xi*{a'^)*\pf\*\Pi\*c'^*(jj'^ — cos{(po)*{ai.pf)*{k.pf)*{k.pi) 

cos{9if)*Ej*Ei*u^—cos{(l)o)*{ai.pf)*{k.pf)*{k.pi)*{a'^)*c'^*u'^+cos{(l)o)*{ai.pf)*{k.pf)*{k.Sf^ 

c^*Ei*uj—cos{(f)Q)*{ai.pf)*{k.pf)*{k.Si)*Xf*Xi*{a'^)*\pf\*c^*Ei*uj'^+cos{(j>o)*{ai.p f)*{k.pf)*{k.Si)*Xf*Xi*{a,'^)*\pi\* 

c'^*cos{9if)*Ef*iv^+cos{(po)*{ai-Pf)*{k-Pi)^*{k.Sf)*Xf*Xi*{a'^)*\pi\*c^*iv — cos{(po)*{ai.pf)* 

Xf*Xi*{a^)*c'^*Ei*(jj—cos{(l)o)*{ai.pi)*{k.pf)'^*{k.pi)*Xf*Xi*\pf\*\pi\*c^—cos{(j)o)*{^^ 

cos{9if)*Ef*Ei-cos{(j)o)*{ai.pi)*{k.pf)^*{k.pi)*c^ + cos{(^a)*{ai-Pi)*{k.pf)'^*{k.^^^ 

{ai.pi)*{k.pfY *{k.piY *Xf*Xi*c^*cos{9if)*Ef*Ei + cos{(f)a)*{ai.pi)*{k.pfY 

{k.pi)*{k.Sf)*Xf*Xi*\pf\*c^*cos{9if)*Ei + cos{(l)i^)*{ai.pi)*{k.pf)'^*{k.pi)*{k.Sf)*Xf*Xi*\pi\*c^ 

{ai.pi)*{k.pf)'^*{k.pi)*Xf*Xi*\pf\'^*c^*coa{9if)*Ei*ij + 2.0*cos{(f)Q)*{ai.pi)*{k.pf)'^*{k.pi)*Xf*Xi*c^ 

Ei*ijj+2.0*cos{(po)*{ai.pi)*{k.pf)'^*{k.pi)*c^*Ef*ij+coa{(po)*{ai.pi)*{k.pf)'^*{k.Si)*Xf*Xi*{a'^)*\pf 

{ai.pi)*{k.pf)*{k.pif*{k.Sf)*Xf*Xi*\pj\*c^*cos{9if)*Ei-cos{(j)Q)*{ai.pi)*{k^^ 

cos{(j)o)*{ai.pi)*{k.pf)*{k.pi)*{k.Si)*Xf*Xi*{a^)*\pf\*c'^*uj + cos{4>o^ 

c'^*w'^—cos{(f)n)*{ai-Pi)*{k.pf)*{k.pi)*Xf*Xi*{a'^)*cos{9if)*Ej*Ei*u'^—cos{(f>n)*{ai.pi)*{k.pf)*{k 

cos{(po)*{ai.pi)*{k.pf)*{k.s f)*{k.Si)*X f*Xi*{a'^)*c^*Ef*oj+cos{(po)*{ai.pi)*{k.pi)*{k.s f)*{k.Si)*Xf*Xi*{a^)*c^*Ef* 

oj+cos{(j)o)*{ai.pi)*{k.pi)*{k.Sf)*Xf*Xi*{a'^)*\pf\*c^*cos{9if)*Ei*Lo'^ — cos{(j)Q)*{ai.pi)* 

\pi\*c'^*Ef*Lo'^ — cos((j)o)*{ai.Sf)*{k.pf)'^*{k.pi)'^*Xf*Xi*\pf\*c^*cos{9if)*Ei — cos{(l)o)*{ai.Sf^ 

Xi*\pi\*c^*Ef+cos{(f)o)*{ai.Sf)*{k.pf)'^*{k.pi)*{k.Si)*Xf*Xi*\pj]*\pi\*c^*cos{9i 

{k.Si)*Xf*Xi*c^^^+cos{(j)o)*{ai.Sf)*{k.pf)'^*{k.pi)*{k.Si)*Xf*Xi*c^*Ef*Ei+cos{(po)*{ai.Sf)*{k.pf^^ 

{a'^)*\Pi\*c'^*(jj—cos{(j)())*{ai.Sf)*{k.pf)'^*{k.Si)*Xf*Xi*{a'^)*c^*Ei*iv+cos{(j)o)*{ai.Sf)*^^ 

cos{9if)*Ei+cos{(i)Q)*{ai.Sf)*{k.pf)*{k.pi)^*Xf*Xi*\pi\*c^*Ef — cos{(po)*{ai.Sf)*{k.pf)*{k.pi)'^*{k.Si)* 

\pi\*c^*cos{9if) — cos{(f)())*{ai.Sf)*{k.pf)*{k.pi)'^*{k.Si)*Xf*Xi*c^" — cos{(f)o)*{ai-Sf)*{k.pf)* 

c^*Ef*Ei — cos{(f)n)*{a-i.Sf)*{k.pf)*{k.pi)'^*Xf*Xi*{a'^)*\pi\*c'^*Lj + 2.0*cos{(j)o)*{ai.Sf)*{k.pf)*{k 

\pi\'^*c^*cos{9if)*ui—2.0*cos{(f)())*{ai.Sf)*{k.pf)*{k.pi)'^*Xf*Xi*\pf\*c'^*cos{9if)*Ef*uj+2.0*cos{(p^ 

{k.pi)'^*Xf*Xi*\pi\*c^*u+cos{(f>())*{ai.Sf)*{k.pf)*{k.pi)*{k.Si)*Xf*Xi*{a'^)*c'^*Ei*Lj+cos{cf)^^^ 

Xf *Xi*{a'^)*\pf\*c'^*cos{9if)*Ei*i^'^ — cos{4>()) *{ai.s f)*{k.pf)*{k. Pi) *Xf*Xi*{a^)*\pi\* c'^*Ef*ij'^ — cos{(po)*{a 

{k.pf)*{k.Si)*Xf*Xi*{a-^)*\pf\*\pi\*c^*cos{0if)*u'^ — cos{<j)o)*{ai.Sf)*{k.pf)*{k.Si)*Xf*Xi*{a'^)*c^^ 

{ai.Sf)*{k.pf)*{k.Si)*Xf*Xi*{a'^)*c'^*Ef*Ei*oj'^+cos{(j)o)*{ai.Si)*{k.pf)^*{k.pi)*Xf*Xi*\^ 

{k.pf)^*{k.pi) *Xf*Xi*\pi\*c^*cos{9if)*Ef — cos{(f)Q) * (ai.Si) * (k .p f)'^ * {k .pi)"^ * X f * Xi* \pf\ *Ei~cos{(j>a) *(ai.Si)* 

{k.pf)^*{k.pi)'^*Xf*Xi*\pi\*c^*cos{9if)*Ef — cos{(po)*{ai.Si)*{k.pf)'^*{k.pi)*{k.Sf)*Xf*Xi*\pf\*\pi\*c^*cos{9if) — 

cos{(j)o)*{ai.Si)*{k.pf)^*{k.pi)*{k.Sf)*Xf*Xi*c^^—cos{(j)o)*{ai.Si)*{k.pf)^*{k.pi)*{k.Sf)*Xf*Xi*c^*Ef*Ei—coa{^^ 

{ai.Si)*(k.pf)'^*{k.pi)*Xf*Xi*{a'^)*\pf\*c^*u!+2.0*cos{(f)o)*(ai.Si)*{k.pf)'^*{k.pi)*Xf*Xi*\pf\'^*\pi\*c^*cos^^ 

2.0*cos{(po)*{ai.Si)*{k.pf)'^*{k.pi)*Xf*Xi*\pf\*c^*oj—2.0*cos{(po)*{ai.Si)*{k.pf)'^*{k.pi)*Xf*Xi*\pi\*c'^*c^^ 

i^+cos{<po)*{ai-Si)*{k.pf)*{k.pi)^*{k.Sf)*Xf*Xi*\pf\*\pi\*c^*cos{6if)+cos{(l)Q)*{ai.Si^ 

c^'^+cos{(f>o)*iai.Si)*(k.pf)*{k.pi)'^*{k.s f)*Xf*Xi*c^*Ef*Ei+cos{(f)Q)*(ai.Si)*{k.pf)*{k.pi)'^*Xf*Xi*{a'^)*\pf\* 
cos{(po)*{ai.Si)*{k.pf)*{k.pi)*{k.Sf)*Xf*Xi*{a'^)*c'^*Ef*(j—cos{(f>o)*{ai.Si)*{k.pf)*{k.pi)*Xf*Xi*{a'^)*\pf 
i^^+cos{(po)* {ai -Si)* {k.pf)*{k.pi)*Xf*Xi*{a'^)*\pi\*c'^*cos{9if)*Ef*Lj^—cos{(f>o)*{ai. Si)* {k.pi)^ 
c'^*Ef*ijj—cos{(f)o)*{ai.Si)*{k.pi)*{k.Sf)*Xf*Xi*{a'^)*\pf\*\pi\*c'^*coa{9if)*iv'^—cos{<po)*{ai.Si)*{k.pi)*{k.Sf)*Xf*Xi* 
{a^)*c^*(^'^+cos{(po)*{ai.Si)*(k.pi)*{k.Sf)*Xf*Xi*(a'^)*c^*Ef*Ei*oj'^—cos{(po)*{k.pf)'^*{k.pi)*Xf*\pf\'^*c'^*cos{^ 
iv*sm(9f)*\a\+cos{(po)*(k.pf)'^*(k.pi)*Xf*\pf\*\pi\*c^*cos{Lpf)*Lu*sm(9i)*\a\ — cos{(po) 
coa{ipf)*cos{9i)*Ef*uj*iiin{9f)*\a\+cos{(j)o)*{k.pf)'^*{k.pi)*Xf*\pi\*c^*cos{Lpf)*cos{9f)*Ef*ij*siu{9i)*\a\ + 
{k.pf)^*{k.pi)*Xf*c^*coii{(fj)*Ej*ijj*sm{9f)*\a\+cos{(j)o)*{k.pf)'^*{k.pi)*Xf*c^*cos{(pf)*Ef*Ei*tJ*sm{9f^ 
cos{(j)o)*{k.pf)'^*{k.pi)*Xi*\pf\*\pi\*c'^*cos{(pf)*uj*sm{9f)*\a\ — cos{(j)o)*{k.pf)^*{k.pi)*Xi*\pf\*c^*cos{(pf)*cos{^ 
Ei*ui*siu{9f)*\a\+cos{(j)o)*{k.pf)'^*{k.pi)*Xi*\pf\*c^*cos{(pf)*cos{9f)*Ei*ij*sm{9i)*\a\ — cos{(po)*{k.pf)'^*{k.p^^ 
Xi*\pi\'^*c'^*cos{(pf)*oj*sm{9i)*\a\+cos{(j)o)*{k.pf)'^*{k.pi)*Xi*c^*cos{(pf)*Ef*Ei*w*sm{9i)*\a\+cos{4^^^ 
{k.pi)*Xi*c^ *cos{(p f)*Ef*ij*sm{9i)*\a\+cos{(po)*{k.pf)*(k.pi)'^*Xf*\pf\'^*c''^*cos{(ff)*ij*sm(9f)*\a]—cos{d^^^ j)* 
{k.pi)'^*Xf*\pf\*\pi\*c'^*cos{(ff)*w*sin{9i)*\a\—cos{(po)*{k.pf)*{k.pi)^*Xf*\pi\*c^*cos{(pf)*cos{9i)*Ef*(^*sm^ 
\a\+cos{(po)*{k.pf)*{k.pi)^*Xf*\pi\*c^*cos{(pf)*cos{6f)*Ef*tJ*sm{6i)*\a\ — cos{<po)*{k.pf)*{k.pi)^*Xf*c^*cos{ipf^ 
Ej*uj*sm{9f)*\a\—coa{(j)Q)*{k.pf)*{k.pi)'^*Xf*c^*coa{(pf)*Ef*Ei*ij*sm{9f)*\a\~cos{<po)*{k.pf)*{k.^^ 

*cos{(pf)*Lv*sm{9f)*\a\—cos{(j)o)*{k.pf)*{k.pi)'^*Xi*\pf\*c^*cos{(ff)*cos{9i)*Ei*(jj*sm{9f)*\a\+cos{^^^ 
{k.pi)^*Xi*\pf\*c^*cos{ipf)*cos{6f)*Ei*u*sm{6i)*\a\+cos{^o)*{k-Pf)*{k.pi)'^*Xi*\pi\'^*c'^*cos{(fif)*w*sm^ 
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cos(0o)*(fc.p/)*(fc.pi)^*Ai*c^*cos((/3/)*£'/*Si*w*sin(6'i)*|o|--cos((/)o)*(fc 

\a\ — cos{(j)Q)*{k.pf)*{k.pi)*Xf*{a'^)*c^*cos{(pf)*Ef*co'^*sin{9f)*\a\ — cos{(j)Q)*{k.pf)*{k.pi)*Xi*{^^ 

w^*sin(6'i)*|a| — cos((/)o)*(A:.p/)*(/s.Si)*Ai*(a^)*|p/|*c^*cos((y9/)*a;^*sin(6'/)*|a|+cos(^o)* 

c^*cos{(ff)*u^*sm{9i)*\a\+cos{(po)*{k.pi)'^*Xf*{a'^)*c^*cos{(ff)*Ef*u^*sm{9f)*\a\+cos{(j)o)*{k 

cos(<^/)*Ej*a;^*sin(^j)*|a|+cos((/)o)*(fc.pi)*(fc.Si)*Ai*(a^)*|p/|*c^*cos((^/)*a;^* 

Aj*(a^)*|pj|*c^*cos((/9/)*w^*sin(6'i)*|a|— cos(0o)*(A;.pi)*A/*(a^)*|p/|^*c^*cos^ 

Xf*{a'^)*\pf\*\pi\*c^*cos{(pf)*u}^*sm{9i)*\a\+cos{(po)*{k.pi)*Xf*{a'^)*\pi\*c'^*cos{(pf)*cos{9i)*Ef* 

cos{(j)o)*{k.pi)*Xf*{a'^)*\pi\*c'^*cos{(pf)*cos{9f)*Ef*u^*sm{9i)*\a\+cos{4>o)*{k.pi)*Xf* 

sin{9f)*\a\—cos{(j)o)*{k.pi)*Xj*{a'^)*c*cos{ipf)*Ef*Ei*uj'^*iiin{9f)*\a\—cos{(j)o)*{^ 
u^*sin{9f)*\a\+cos{(po)*{k.pi)*Xi*{a'^)*\pf\*c'^*cos{ipf)*cA)s{9i)*Ei*ij^*s\n{9f^^ 
c'^*cos{(pf)*cos{9f)*Ei*u^*sm{9i)*\a\+cos{4>o)*(k:pi)*Xi*{a'^)*\pi\'^*c^*cos{^ 
Xi*{a'^)*c*cos{(ff)*Ef*Ei*Lu^*sm{9i)*\a\—cos{(pQ)*{k.pi)*Xi*{a^)*c*cos{(ff)*Ef*uj^*sin{9i)*\a^ 
{k.pf)^*{k.pi)^*Xf*Xi*\pf\*\pi\*c^+sm{(j)o)*{a2-Pf)*{k.pf)^*{k.pi)^*Xf*Xi*c^*cos{9if)*Ef*Ei+sm^ 
{k.pf)'^*{k.pi)'^*c^—sm{(j)o)*{a2-Pf)*ik.pf)'^*{k.pi)*{k.Si)*Xf*Xi*\pf\*c^*Ei—sm{(j)^^ 
Xf*Xi*\pi\*c'^*cos{9if)*Ef-sm{(po)*{a2.pf)*{k.pf)*{k.pi)^*Xf*Xi*\pf\*\pi\*c^-sm{^^^ 
Xi*c'^*cos{9if)*Ef*Ei—sm{(po)*{a2-Pf)*{k.pf)*{k.pi)^*c^+sin{(j)o)*{a2.pf)*{k.pf)*{k.pi)'^*{^ 
sm{(j)o)*{a2-Pf)*{k.pf)*{k.pi)^*{k.Si)*Xf*Xi*\pi\*c^*cos{9if)*Ef — 2.0*sm{4io)*{a2.pf)*{k.pf)*{k.p 
c^*cos{9if)*Ef*iv + 2.Q*sm{(t)o)*{a2-Pf)*{k.pf)*{k.pi)'^*Xf*Xi*c'^*cos{9if)*Ef*Ef*ij + 2.Q*sm^ 
{k.pi)'^*c^*Ei*iv — sm((f>Q)*{a2.pf)*(k.pf)*{k.pi)*(k.Sf)*Xf*Xi*{a'^)*\pi\*c^*iv+si^^^ 
Xi*{a^)*\pf\*\pi\*c^*ij'^ — iim{4>Q)*(a2-Pf)*{k.pf)*{k.pi)*Xf*Xi*{a'^)*coii{9if)*Ef*Ei* 

{k.pi)*{a'^)*c'^*LV^+sm{(j)o)*{a2-Pf)*{k.pf)*{k.s f)*{k.Si)*Xf*Xi*{a'^)*c'^*Ei*ui — iiin{(t>o)*{a2-Pf)*{k.pf)*{k.Si)*Xf* 

Xi*{a'^)*\pf\*c'^*Ei*u^+sin{4>o)*(a2-Pf)*(k.pf)*(k.Si)*Xf*Xi*(a^)*\pi\*c^*cos^^ 

{k.Sf)*Xf*Xi*{a^)*\pi\*c^*ij—sm{(j)o)*(a2-Pf)*(k.pi)*(k.Sf)*{k.Si)*Xf*Xi*{a.'^ 

{k.pi)*Xf*Xi*\pf\*\pi\*c^—sm{(j)Q)*{a2.pi)*{k.pf)'^*{k.pi)*Xf*Xi*c'^*cos{9if)*Ef*Ei 

c^+sin(^o)*(o2-Pi)*(fc-P/)^*(fc-Pi)^*A/*Ai*|p/|*|pj|*c^+sin((^o)*(«2-Pi)*(fc._p/)^*(fc 

Ei+sin{<po)*{a2-Pi)*{k.pf)'^*{k.pi)'^*c^+sm((po)*{a2.pi)*{k.pf)'^*{k.pi)* 

{a2-Pi)*{k.pf)^*{k.pi)*{k.Sf)*Xf*Xi*\pi\*c^*Ef — 2.0*sm{(j)o)*{a2.pi)*{k.pf)^*{k.pi)*Xf*Xi*\pf\'^*c^ 
u+2.0*sm{^o)*{a2.Pi)*{k.pf)'^*{k.pi)*Xf*Xi*c'^*cos{6if)*Ej*Ei*u+2.0*sm{(j)o)* 

sm{(j>o)*ia2-Pi)*(k.pf)'^*(k.Si)*Xf*Xi*{a'^)*\pf\*c'^*iv—sm{(t>o)*ia2-Pi)*{k.pf^^ f)*Xf*Xi*\pf\*c^*cos{9if)* 

Ei-sm{<j)o)*{a2.Pi)*{k.pf)*(k.pi)'^*{k.Sf)*Xf*Xi*\pi\*c^*Ef-sm{(j>o)*{a2-P^ 

c'^*u+sm{(j)o)*{a2.pi)*{k.pf)*{k.pi)*Xf*Xi*{a'^)*\pf\*\pi\*c'^*uj'^—sm{(l)o)*{a2-^^^ 

Ef*Ei*oj'^ — sm((f>Q)*{a2.pi)*{k.pf)*{k.pi)*{a'^)*c**U!'^ — sm(i^Q)*{a2.pi)*{^^ 

uj+sm{(j)o)*{a2-Pi)*(k.pi)*{k.Sf)*{k.Si)*Xf*Xi*{a^)*c**Ef*Li!+iiin((^Q)*{a2.pi^ 

cos{9if)*Ei*u'^ — sin{4>o)*{a2.pi)*{k.pi)*{k.Sf)*Xf*Xi*{a'^)*\pi\*c'^*Ef*ij^—sin{(po)*{a^ 

\pf\*c^*cos{9if)*Ei—sm{(j)o)*{a2.Sf)*{k.pf)^*{k.pi)^*Xf*Xi*\pi\*c^*Ef+sm{(j)o)*{a2-Sf)*{k.pf^^ 

Xi*\pf\*\pi\*c^*cos{9if) + sm{(pa)*{a2-Sf)*{k.pf)'^*{k.pi)*{k.Si)*Xf*Xi*c^^+sin{^^ 

Xf*Xi*c^*Ef*Ei+sin{(j)o)*{a2.Sf)*{k.pf)'^*{k.pi)*Xf*Xi*{a^)*\pi\*c^*Lv—iiin{^^^ 

c *i?j*w+sin(0o)*(a2-s/)*(A;.p/)*(A;.pip*A/*Ai*|p/|*c^*cos(0i/)*i?i+s 

Ef—sm{(po)*{a2-Sf)*{k.pf)*{k.pi)'^*{k.Si)*Xf*Xi*\pf\*\pi\*c^*cos{Oif)—sm{(f)o)*{^^ 

c^^—sin{(f>o)*{a2-Sf)*{k.pf)*{k.pi)'^*(k.Si)*Xf*Xi*c^*Ef*Ei—sm{(j)o)*{a2-Sf)*{k.pf)* 

2.0*sm{(j)o)*{a2-Sf)*{k.pf)*{k.pi)^*Xf*Xi*\pf\*\pi\'^*c^*cos{9if)*ij—2.0*sin{(j)o)*{a2.Sf)^ 

c'^*cos{6if)*Ef*u+2.0*sm{(t)Q)*{a2.Sf)*{k.pf)*{k.pi)^*Xf*Xi*\pi\*c^*ivsm{(j)o)*{^^ 

{a'^)*c'^*Ei*uj+sm{<po)*{a2-s f)*{k.pf)*{k.pi)*Xf*Xi*{a^)*\pf\*c^*cos{9if)*Ei*oj'^—sm^ 

Xi*{a'^)*\pi\*c^*Ef*Lu'^—sm{(t>o)*(a2-Sf)*{k.pf)*(k.Si)*Xf*Xi*{a^)*\pf\*\pi\*c**cos{^ 

{k.Si)*Xf*Xi*{a'^)*c^*iv'^+sin((f>o)*{a2.Sf)*(k.pf)*(k.Si)*Xf*Xi*{a^)*c^*Ef*Ei*Lu'^ 

Xf*Xi*\pf\*c^*Ei+sm{(j)o)*{a2-Si)*{k.pf)^*{k.pi)*Xf*Xi*\pi\*c^*cos{9if)*Ef—sin{(j)o)*{a2-Si^ 

Ai*|p/|*c''*i?i-sin((/)o)*(a2.Si)*(/c.p/)2*(fc.pj)2*A/*Ai*|pi|*c''*cos(%)*i?/-sin 

Xf*Xi*\pf\*\pt\*c^*coii{9if)-sm((f>o)*{a2.Si)*(k.pf)'^*(k.pi)*(k.Sf)*Xf* 

Xf*Xi*c^*Ef*Ei—sm{(j)o)*{a2.Si)*{k.pf)'^*{k.pi)*Xf*Xi*{a'^)*\pf\*c^*w+2.0*sm{(po^ 

\pf\^*\pi\*c^*cos{Oif)*tJ+2.0*sm{(l)o)*{a2.Si)*{k.pf)^*{k.pi)*Xf*Xi*\pf\*c^*tJ—2.0*s^^ 

Xf*Xi*\pi\*c'^*cos{9if)*E'j*Lu+sm(4)())*{a:2--'^i)*ik.pf)*{k.pi)'^*{k.Sf)*Xf* 

{k.pf)*{k.pi)^*{k.Sf)*Xf*Xi*c^^sm{(j)o)*{a2.Si)*{k.pf)*{k.pi)^*{k.Sf)*Xf*Xi*c^*Ef*Ei + sin{^^^^ 
{k.pi)^*Xf*Xi*{a^)*\pf\*d^*i^+sm{(po)*{a2-Si)*{k.pf)*{k.pi)*{k.Sf)*Xf*Xi*{a'^)*d^*E 

{k.pi)*Xf*Xi*{a'^) * \pf\*c'^*Ei*uj'^ + s'm{(f>Q) * {a2.Si) *{k.pf)* {k.pi)*Xf*Xi * (a^) * \pi\*c'^*cos{9if)*Ef*uj'^—sm{(f>Q) * 

{a2-Si)*{k.pi)^*{k.Sf)*Xf*Xi*{a'^)*c'^*Ef*ij~sm{(po)*{a2-Si)*{k.pi)*{k.Sf)*Xf*Xi*{a'^)*\pf\*\p^^ 

sm{(j)o)*{a2-Si)*{k.pi)*{k.Sf)*Xf*Xi*{a'^)*c^*uj'^+sm{(j)())*{a2-Si)*{k.pi)*{k.Sf)*Xf*Xi*{a'^)* 

{k.pf)'^*{k.pi)*Xf*\pf\'^*c''*ij*sin{(pf)*sm{9f)*\a\+sin((f)Q)*{k.pf)'^*{k^^^ 

sm{(j)o)*{k.pf)^*{k.pi)*Xf*\pi\*c^*cos{9f)*Ef*uj*sm{ipi)*sm{9i)*\a\—sm{(^^^^ 
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uj*sm{(pf)*sm{0f)*\a\+sin((f)o)*{k.pf)'^*{k.pi)*Xf*c^*E'j-*iv*sm{^ 

Ei*w*sm{(fif)*sm{9f)*\a\+sin{(po)*{k.pf)'^*{k.pi)*Xi*\pf\*\pi\*c^*u*sm{ipf)*si^^ 

Xi*\pf\*c^*cos{6f)*Ei*u*sin{ipi)*sm{Oi)*\a\—sm{(po)*{k.pf)'^*{k.pi)*Xi^^ 

sm{(j)o)*{k.pf)'^ *{k.pi)*Xi*\pi\'^*c''' *iv*sm{ifi)*sm{9i)*\a\+sm{^ 

\a\+sin{(l)o)*{k.pf)^*{k.pi)*Xi*c^*Ef*iv*sin{ifi)*sm{9i)*\a\+sm{(j)o)*{k.pf^ 

\a\—sm{(j)o)*{k.pf)*{k.pi)'^*Xf*\pf\*\pi\*c'^*uj*sm{ipi)*sin{9i)*\a\+sin{(po)*{k.pf^ 

uj*sm{(pi)*sm{6i)*\a\—sm{(j)o)*(k.pf)*{k.pi)'^*Xf*\pi\*c^*cos{9i)*Ef*Lu*sm{c^ 

Xf*c^*E'j*(M*sin{(pf)*sm{9f)*\a\^sm((}>Q)*{k.pf)*{k.pi)'^*Xf*c^*Ef*Ei*uj*sm{(ff)* 

{k.pi)'^*Xi*\pf\*\pi\*c'^ *uj*sm{ipf)*sm{9f)*\a\+sm{(l)o)*{k.pf)*{k.pi)'^*Xi*\pf\ 

\a\—sm{(j)o)*(k.pf)*{k.pi)'^*Xi*\pf\*c^*cos{9i)*Ei*iv*sm{(pf)*sm{9f)*\a\+sm{(^^ 

sm{(pi)*sm(9i)*\a\—sin((pQ)*{k.pf)*{k.pi)'^*Xi*c^*Ef*Ei*oj*sm{(pi)*s'^^^^ 

u}*sin{ipi)*sm{9i)*\a\—sin{(po)*{k.pf)*{k.pi)*Xf*{a'^)*c^*Ef*i^'^*sin((pf)*^^ 

C^*Ei*ij^*sin{(fii)*sm{9i)*\a\—sin{(po)*{k.pf)*{k.Si)*Xi*{a'^)*\pf\*c^*uj'^*sin^ 

Xi*{a'^)*\pi\*c^*u^*sm{ipi)*sm{9i)*\a\+sm{(po)*{k.pi)'^*Xf*(a'^)*c^*Ef*ijj^ 

{a^)*c^*Ei*u'^*sm{Lpi)*sm(9i)*\a\+sin{(po)*(k.pi)*(k.Si)*Xi*(a'^)*\p 

{k.Si)*Xi*{a'^)*\Pi\*c^*ui'^*sm{(pi)*sm{9i)*\a\—sm{(j)o)*{k.pi)*Xf*{a'^)*\pf\'^*c^ 

{k.pi)*Xf*{a^)*\pf\*\pi\*c^*cu^*sm{(fi)*sm{9i)*\a\—sm{(po)*{k.pi)*Xf*{a^ 

\a\+sm{4>o)*{k.pi)*Xf*{a'^)*\pi\*c'^*cos{9i)*Ef*iv^*sin{(ff)*sm(9f)*\a\+sm{(^^ 

sm{6f)*\a\ — sm{(/)o)*{k.pi)*Xf*{a'^)*c*Ef*Ei*u^*sm{(fif)*sm{9f)*\a\ — sm{(j)o)*{k 
sm{(ff)*sin{9f)*\a\—sm{(j)o)*{k.pi)*Xi*{a'^)*\pf\*c'^*cos{9f)*Ei*u^*sm^ 

*cos{9i)* Ei*ij'^ *sm(Lp f) *siii{9 f)*\a\+sin{(f)o)*{k. Pi)* Xi*(a'^)*\pi^ 

{a^)*c*Ef*Ei*u^*sin{ipi)*sin{9i)*\a\—sm{(j)o)*{k.pi)*Xi*{a'^)*c*Ef*Lo^*s^ >/[2.0*{k.pf)'^*{k.pi)'^*c'^] 
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